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Abstract 

We study a system of two pointlike particles coupled to three dimensional Einstein grav- 
ity. The reduced phase space can be considered as a deformed version of the phase space 
of two special-relativistic point particles in the centre of mass frame. When the system is 
quantized, we find some possibly general effects of quantum gravity, such as a minimal dis- 
tances and a foaminess of the spacetime at the order of the Planck length. We also obtain a 
quantization of geometry, which restricts the possible asymptotic geometries of the universe. 



Outline and summary 



The Kepler system is the simplest realistic example of a coupled two body system, and belongs to 
the few systems that can be solved exactly within the framework of Newtonian gravity. It consists 
of two pointlike objects, characterized only by their masses, and interacting with the gravitational 
field. Unfortunately, within the framework of general relativity, the two body problem not only 
lacks of an exact solution. It is not even well defined, because Einstein gravity in four spacetime 
dimensions does not admit pointlike matter sources. Clearly, this makes general relativity so 
interesting. But the obvious drawback is the absence of a simple but still realistic toy model, 
which is sometimes very useful. The Kepler system is, in a sense, the hydrogen atom of gravity. 

The situation is different in three spacetime dimensions, where Einstein gravity is not only a 
much simpler field theory [|], ^ [3j]. It also admits pointlike matter sources 0, ^ 0, H ^ |lOl ], 
and even a more or less straightforward canonical quantization [II], 0, [TTl [T2|, O). The vacuum 
Einstein equations in three dimensions require the spacetime to be flat outside the matter sources. 
There are neither gravitational waves, nor local gravitational forces. However, the spacetime 
becomes curved if matter is present. The simplest example of a non-trivial spacetime is the 
gravitational field of a massive pointlike particle. The spacetime is the direct product of a real 
line with a conical space. The particle is sitting at the tip, and the deficit angle of the cone is 
87rGm, where m is the rest mass of the particle, and G is Newton's constant. 

In units where the velocity of light is one, G has the dimension of an inverse mass or energy. 
The inverse M P1 = 1/G is the Planck mass. It is a classical quantity in the sense that % is not 
involved in the definition. There is thus a certain amount of curvature on the world line of a point 
particle, which is proportional to its mass. But this is just a very simple conical singularity in 
an otherwise flat spacetime. There is a simple way to visualize a spacetime containing a single 
point particle. One starts from a flat three dimensional Minkowski space and cuts out a wedge. 
A wedge is a subset which is bounded by two timelike half planes, whose common boundary is 
a timelike geodesic. This geodesic becomes the world line of the particle. 

The two half planes are mapped onto each other by a certain isometry of Minkowski space. 
It is a Lorentz rotation about the world line, and the angle of rotation is RnGm. The region 
inside the wedge is taken away, and the points on the two half planes are identified, according 
to the Lorentz rotation. The result is a locally flat spacetime with a conical singularity on the 
world line. For a massless particle, the same procedure can be applied to a pair of half planes, 
whose common boundary is a lightlike geodesic. They are then mapped onto each other by a 
null rotation. The Kepler spacetime contains two particles, and therefore we have to apply this 
procedure twice. The result is shown in figure |T} 

Unless one of the particles has a deficit angle which is bigger than n, hence a mass above 
Mpi/8, every possible two particle spacetime can be constructed in this way. The geometry of the 
spacetime only depends on the relative motion of the particles, which can be read off immediately 
from the relative orientation of the two world lines in Minkowski space. It is therefore possible to 
get an overview of all possible spacetime geometries very easily. However, there are also some 
problems with this simple construction. First of all, it does not work for masses bigger then 
Mpi/8, although such spacetimes do exist. In this case, the half planes defining the boundary of 
the wedges must be replaced by curved surfaces, as otherwise the wedges overlap. But this is 
actually not a serious problem. 
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Figure 1: The Kepler spacetime can be constructed by cutting out two wedges from a fiat 
Minkowski space. The faces are identified, such that two conical singularities arise in an 
otherwise flat spacetime. In the rest frame of the each particle, the deficit angle of the conical 
space is proportional to the mass of the particle. 



The more serious problem has to do with the asymptotic structure of the spacetime at infinity. 
The region far away from the particles is split into two segments in figure [j]. Each segment is a 
subset of Minkowski space. But on the wedges we have to apply non-trivial transition functions, 
relating the Minkowski coordinates on one side to those on the other side. To find out what the 
spacetime looks like at infinity, it would be nicer to have a single coordinate chart covering this 
region. There is in fact a particular reason why we are interested in the asymptotic structure 
of the Kepler spacetime. In order to quantize it in the end, we first have to set up a proper 
classical Hamiltonian formulation. This requires a proper definition of an action principle for the 
underlying field theory of Einstein gravity. And this again requires some kind of asymptotical 
flatness condition to be imposed on the metric at infinity fll4{]. 

The asymptotic structure of the Kepler spacetime depends crucially on the relative motion of 
the particles. If they are moving slowly, then far away from the particles the spacetime is also 
conical. It looks almost like the gravitational field of a single particle, whose mass is equal to 
the sum of the two masses of the real particles. The rest frame of this fictitious particle can be 
identified with the centre of mass frame of the universe. If the particles are moving faster, the 
apparent mass of the fictitious particle has to be replaced by the total energy of the system. It 
also receives a spin, which represents the total angular momentum. But still, the universe looks 
like a cone at infinity, and this cone defines the centre of mass frame. 

Something strange happens when the relative motion of the particles exceeds a certain thresh- 
old [|T5|]. The definition of a centre of mass frame then breaks down, and the asymptotic structure 
of the spacetime is no longer conical. Even more peculiar, the spacetime then contains closed 
timelike curves [ [T6| , [FT], [T8|]. Clearly, these are very interesting features of such a simple two 
particle spacetime. But for our purpose we have to exclude them, again because we want to set 
up a proper Hamiltonian framework. This requires a well defined causal structure of the space- 
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time. Otherwise the Hamiltonian, or ADM formulation cannot be applied to general relativity. It 
requires the existence of a foliation of the spacetime by spacelike slices [ JT9| ] 

In order to get rid of these problems, we have to go over from the simple construction of 
the Kepler spacetime in figure Q], to a slightly more sophisticated one, which focuses more on 
the asymptotic structure of spacetime at infinity. The basic idea is, first to fix the asymptotic 
structure of the spacetime, and then insert the particles. In figure [IJ the sequence is the other 
way around. The world lines are inserted first, then the actual spacetime is constructed, and 
finally the asymptotic structure can be read off by looking at the region at infinity. The transition 
from this picture to an alternative description of the Kepler spacetime, based on its asymptotic 
structure, is explicitly carried out in p^ ] . There we also give a comprehensive overview of all 
those spacetimes that admit the definition of a centre of mass frame. 

Somewhat schematically, the alternative construction of a two particle spacetime is shown in 
figure One starts from a big cone, cuts off the tip, and identifies the cut lines, which are 
two geodesies, such that a conical surface with two tips arises. The three dimensional version 
of this construction yields the Kepler spacetime, although the technical details are a little bit 
more involved. The advantage of this procedure is that the original cone immediately defines the 
centre of mass frame of the universe, as seen by an observer at infinity, and independent of the 
way the particles are inserted. Moreover, we will be able also introduce position and momentum 
coordinates of the particles, referring to the centre of mass frame defined by the big cone. 

In this way, the Kepler system can effectively be treated like a simple two particle system in a 
fixed three dimensional background spacetime, although the technical details are again slightly 
more involved. The phase space becomes a finite dimensional manifold, and the kinematical and 
dynamical properties of the Kepler system are finally encoded in the usual way, in the symplectic 
structure, or the Poisson bracket, and the Hamiltonian. The phase space structure is very similar 
to a system of two special-relativistic point particles in a flat Minkowski space, restricted to the 
centre of mass frame. And in fact, it is possible to take the limit G — > 0, where the gravitational 
interaction is switched off, and then the Kepler system reduces to this free particle system. 

The article is therefore organized as follows. In section [1], we shall only study the free particle 
system. Of course, there is nothing new about this to be learned, expect perhaps the relativistic 
definition of a centre of mass frame, which is in this form not a very common concept. It is 
actually a centre of energy frame, but we shall stick to the more familiar notion of a centre of 
mass. After imposing the appropriate restriction on the phase space, we shall go through the usual 
canonical programme. Starting from the classical Hamiltonian framework, we perform a phase 
space reduction, derive and solve the classical equations of motion, and finally we quantize it, 
deriving the energy eigenstates and the spectra of certain interesting operators. The idea behind 
this preparation, using a well known and very simple toy model for our real toy model, is to keep 
the conceptual aspects apart from the technical aspects. 

The conceptual aspects are, for example, the definitions of the various phase spaces, the way 
the mass shell constraints are imposed, the principle idea of the phase space reduction, and fi- 
nally also the quantization methods. At the classical level, the free particle system also provides a 
nice toy model for the Hamiltonian formulation of general relativity within the ADM framework. 
This is indicated in figure 0. At the quantum level, we are going to consider two alternative quan- 
tization methods, the Schrodinger method applied to the reduced classical phase space, where all 
gauge symmetries are removed, and the Dirac method applied to an extended phase space, where 
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the dynamics of the system is defined by a generalized mass shell constraint. 

All these concepts can then be applied to the Kepler system in the very same way. It is there- 
fore useful first to explain them using a much simpler model. The more technical aspects are 
then the modifications that we have to make when the gravitational interaction is switched on. 
They are divided into a classical part in section 0, and a quantum part in section |3|. In the clas- 
sical section, we shall first look at the Kepler spacetime itself, in the way explained above and 
described in figure |3|. We shall mainly focus on a geometric point of view, without going into 
any technical details and proofs. All these details can be found in the previously mentioned and 
more comprehensive article p0[]. It includes, in particular, the precise transition from figure [I] to 
the alternative description in figure [5[ 

In a certain sense, the Kepler spacetime is a deformation of the free particle spacetime, with 
Newton's constant being the deformation parameter. At any stage of the derivation we can always 
get back to the free particle system by taking the limit G — > 0. This can be used as a cross check 
at various points, and sometimes it is even possible to set up general rules, telling us how to 
deform the free particle system to obtain the corresponding structures of the coupled system. 
This also applies to the phase space of the Kepler system, which is the subject of the second 
part of section It is a deformed version of the free particle phase space. The interesting 
point is thereby that gravity affects the symplectic structure of this phase space, and not only the 
Hamiltonian, which is otherwise the typical feature of interactions between point particles. 

The actual derivation of the deformed symplectic structure will not be given in the article. 
We should however emphasize that this derivation is a very important point. We do not want 
to make any special or unmotivated assumptions. Instead, the only assumption that enters the 
definition of the Kepler system is the following. The kinematical and dynamical features of the 
gravitational field are completely defined by the Einstein Hilbert action. All the relevant phase 
space structures can then be derived from this action principle, by a straightforward phase space 
reduction. However, apart from the symplectic structure, all other features of the phase space 
can more or less be inferred from geometric considerations. We shall therefore restrict to these 
geometrical aspects here, and refer to [ ^0| ] for the derivation of the symplectic structure for a 
more general multi particle system. 

In the last part of section Q, we will go through the whole canonical programme once again, 
deriving and solving the classical equations of motion, and briefly describing the various kinds of 
trajectories. At this point, we can actually forget about the general relativistic nature of the sys- 
tem, and treat it as if it was a simple two particle system living in a three dimensional background 
spacetime. Or, if we do not want to give up the general relativistic point of view completely, we 
may at least stick to the ADM picture, and consider the Kepler system as a space which evolves 
in time. As shown in figure [|, the phase space variables define the geometry of space at a mo- 
ment of time, and this geometry changes with time. Effectively, the particles are moving in a two 
dimensional space. We can also think of a typical scattering process, and define quantities like 
incoming and outgoing momenta, and scattering angles. 

In the quantum section, we will first try to apply the same quantization methods that we also 
applied to the free particle system. We'll find that the straightforward Schrodinger quantization 
fails, due to some peculiar features of the deformed classical phase space and its symplectic 
structure. This already indicates that there are some new effects to be expected, which are due to 
the gravitational interaction, and which are fundamentally different from other interactions. The 
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Dirac method however works. It is possible to set up a well defined operator representation, and 
to quantize and solve the constraint equation, which is a generalized Klein Gordon equation. A 
similar equation has also been found for a somewhat simpler single particle system [|T3|]. 

We can solve this constraint equation, and finally we are able to express the energy eigenstates 
of the Kepler system explicitly as wave functions on a suitably defined configuration space. Suit- 
ably thereby means that the wave function has the usual physical interpretation as a probability 
amplitude for the particles in space. Or, once again, if we want to stick to the general relativistic 
point of view, it is a probability amplitude for certain geometries of space. We have in this sense 
a simple example for a truly quantized general relativistic system. The quantum state gives us 
probabilities for geometries. But nevertheless, it is somewhat more intuitive to consider the wave 
function as a probability amplitude in an ordinary space. The space is thereby the configuration 
space of the two particles in the centre of mass frame. And it is also useful to have in mind the 
usual picture of a scattering state in quantum mechanics. 

The energy eigenstates are parameterized by two quantum numbers. The inverse radial wave- 
length at infinity represents the eigenvalue of the incoming and outgoing momenta of the parti- 
cles. And the inverse angular wavelength represents the eigenvalue of the angular momentum. 
The latter is quantized in steps of h, and the radial momentum turns out to have a positive con- 
tinuous spectrum. The energy, or the frequency of the wave function, only depends on the radial 
quantum number. So far, these are the typical features of scattering states in quantum mechan- 
ics. However, the energy spectrum is bounded from below and from above. This is actually not 
surprising, because the total energy contained in a three dimensional universe is bounded from 
above by M P1 /4, which corresponds to the maximal deficit angle 2% of a conical spacetime. The 
lower bound for the total energy is the sum of the rest masses of the particles, which is also not 
surprising. 

What is remarkable, however, is that the radial momentum of the particles is nevertheless 
unbounded. For very large momenta of the particles, the total energy of the Kepler system ap- 
proaches the upper bound Mpi/4. This behaviour has also been found for a single particle, where 
the upper bound is half as big, thus Mpi/8 Jl3t]. The relation between the spatial momentum of 
the particles in the centre of mass frame, and the total energy of the system is shown in figure H, 
where it is compared to the corresponding free particle energy. The energy of the Kepler system 
is always smaller than the free particle energy. The difference is a kind of gravitational binding 
energy. At the quantum level, this has the strange consequence that the wavelength of the wave 
function in space can be arbitrarily small, but the frequency is bounded from above. 

Finally, we shall then look at the wave functions themselves. Unfortunately, it is hardly pos- 
sible to read off any physically interesting information directly from the analytic expressions. 
They are somewhat complicated, involving hypergeometric functions. We shall therefore look 
at the graphical representations of some typical wave functions. They are shown in the figures |^ 
through both for the free particles and the Kepler system. We shall thereby find the following 
interesting features. At large distances, the gravitational interaction has almost no effect. We 
have the typical scattering wave function, which is a superposition of an ingoing and an outgoing 
radial wave. The only difference between the free and the coupled system is a phase shift, which 
indicates that some interaction takes place when the particles are closer to each other. 

At small distances however, which are of the order of the Planck length i = Gh, the wave 
function changes drastically. We typically find that the particles avoid to be at the same point 
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in space. We even find that under certain very general circumstances, it is impossible for the 
particles to get closer to each other than a certain minimal distance. It is of the order of ten to 
hundred Planck lengths, and depends on the rest masses of the particles and their statistics. In 
three spacetime dimensions, there are not only bosons and fermions, but also anyons, and there 
is also a generalized statistics if the particles are not identical [|2~T|]. All this can be taken into 
account very easily when the quantization is performed. The finite lower bound for the distance 
of the particles in space arises whenever the particles are not two bosons, and it is maximal for 
two identical fermions. 

Referring again to the scattering picture, it is reasonable to say that due to this minimal dis- 
tance, it is impossible to probe the structure of spacetime at small length scales, even if we 
increase the momentum of the particles unboundedly. This is exactly the kind of limit that quan- 
tum gravity is expected to impose on the ability to look at small length scales in spacetime. 
Another feature of our toy model is closely related to this, but a little bit more general. Even 
if the particles are further apart than the minimal distance, it is still impossible to localize them 
within a box that is smaller than a certain size. More precisely, it is impossible to find a quantum 
state where the relative position of the particles in space is arbitrarily sharp at a given moment 
of time, even if the distance between the particles is many orders of magnitude above the Planck 
scale. 

Both features indicate that the quantized spacetime in which the particles are living obtains a 
kind of foamy structure. Unfortunately, it is not possible to derive a more explicit and intuitive 
spacetime spectrum, like the one for the single particle system in |JT3|] . But in principle, we have 
a very similar situation, and this is also expected to arise in a more realistic, or even in a fully 
consistent theory of quantum gravity in higher dimensions. Hence, although the Kepler system 
is only a very simple toy model, some principle effects of quantum gravity can be seen. These 
are really quantum gravity effects, because they disappear not only at the classical level, hence in 
the limit % — > 0, whatever this precisely means, but also in the well defined limit G — > 0, where 
the gravitational interaction is switched off. 

Finally, a nice feature of this toy model is that, once the reduction to a two particle system 
in three dimensions is accepted, everything else can be derived exactly and without any further 
assumptions. It is possible to keep the assumptions and simplifications clearly apart from the 
mathematics and the physical conclusions. There are no hidden points were additional, say, 
intuitive assumptions must made. The only assumption that enters the definition of the Kepler 
system as a toy model is the Einstein Hilbert action with the appropriate matter terms for the 
particles, which is assumed to define the dynamics of the gravitational field. This is the concept 
on which the derivation of a general multi particle phase space is based in [JTOQ. The point where 
this assumption enters this article, is the definition of the symplectic structure in section ^|. 



1 The free particle system 

Before we switch on the gravitational interaction, let us consider a system of two uncoupled 
relativistic point particles iXk (k = 1,2) in flat, three dimensional Minkowski space. As a vector 
space, we identify this with the spinor representation sl(2) of the three dimensional Lorentz 
algebra. The twelve dimensional kinematical phase space is spanned by the positions x k = x^ja 
and the momentum vectors pk = plla of the particles, where "y a (a = 0, 1, 2) is an orthonormal 
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basis given by the usual gamma matrices ( |A.1| ). Further conventions and some useful formulas 
regarding the vector and matrix notation are given in the appendix. 

We have the usual symplectic potential 0, from which we read off the Poisson brackets, 



6> = i£Tr(p fc d* fc ) =► {p a k ,x b k } = V ab . (1.1) 
k 

The Hamiltonian if is a linear combination of the two mass shell constraints with Lagrange 
multipliers Ck € M as coefficients, 

H = Y,CkC k , C k = \ r Ti{p k 2 ) + \m k 2 . (1.2) 

k 

The physical phase space is a subset of the kinematical phase space, which is defined by the mass 
shell constraints and the positive energy conditions, 

±Tr(p fe 2 ) = -m k \ p\ = |Tr(p fc7 °) > 0. (1.3) 

The world lines are parameterized by a common, unphysical time coordinate t, and the Hamilto- 
nian generates the time evolution with respect to this coordinate, 

p k = {H,p k } = 0, x k = {H,x k } = (kp k . (1.4) 

The freedom to choose the multipliers (k corresponds to the gauge freedom to reparameterize 
the world lines. Finally, there are some conserved charges which are of interest, namely the total 
momentum vector and the total angular momentum vector, 

p = Y.Pki J = \Y}p^ x k}- (i.5) 

k k 

The associated rigid symmetries are the translations and Lorentz rotations of the world lines with 
respect to the reference frame, which is defined by the coordinates of the embedding Minkowski 
space. 

The centre of mass frame 

So far, this is the standard Hamiltonian formulation of a special-relativistic two particle system. 
Since we are only interested in the relative motion of the particles, we shall now impose some 
further restrictions on the phase space variables. Provided that the total momentum vector is 
positive timelike, there always exists a reference frame where 

P = Af 7o , J = 5 7o , M,SeR. (1.6) 

We call this the centre of mass frame. There is only one special situation where a centre of mass 
frame does not exist. If both particles are massless and if they move with the velocity of light 
into the same direction, then the total momentum is lightlike. These special states are excluded 
in the following. 

In the centre of mass frame, the rigid symmetries are reduced to a two dimensional group of 
time translations and spatial rotations about the 70-axis. The associated charges are the total 



7 



Figure 2: The free particle system in the centre of mass frame. The embedding Minkowski 
space is foliated by a family of equal time planes, labeled by an ADM time coordinate t. The 
relative position of the particles in space is defined by X and <fi, and the spatial momentum by 
K and j3. The energies of the particles are and the clock T represents the absolute time 
in the centre of mass frame, which can be regarded as a reference frame of some external 
observer. 



energy M and the spatial angular momentum S. It is also allowed to speak about an absolute 
time, defined by the 7 -axis, and an absolute space orthogonal to it, defined by the 7i )2 -axes in 
Minkowski space. Essentially, we have a non-relativistic system, although the dynamical prop- 
erties of the particles are still relativistic. We can use this to impose a gauge condition, which 
restricts the way the world lines are parameterized. It is reasonable to choose the parameteriza- 
tion such that 

x° l = x° 2 . (1.7) 

At each moment of time t, the particles are then located on the same equal time plane, as indicated 
in figure |2|. We can think of an ADM like foliation of the embedding Minkowski space by equal 
time planes. The planes are labeled by an ADM time coordinate t, and each plane represents an 
instant of time in the centre of mass frame. For the moment, we do not require the ADM time 
t to be related in any way to the absolute time in the centre of mass frame. Hence, there is still 



one gauge degrees of freedom left, which is compactible with the gauge condition (J1.7[). This is 
a simultaneous reparameterization of both world lines. 

From the phase space point of view, the various restrictions can be regarded as additional 
constraints. All together, we have seven constraints. The two mass shell constraints (|1.3[), the 
gauge condition (|T77|), and four independent spatial components of the definition ( |1.6p of the 
center of mass frame. Note that the time components of these equations are definitions of the 
phase space functions M and S. Since we have an odd number of constraints, it is clear that at 
least one of them is a first class constraint, which generates the simultaneous time evolution of 
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the two particles as a gauge symmetry. The other six are second class constraints, as we are now 
going to show. 

To easiest way to do this is to eliminate six of the seven constraints, and to show that the 
reduced symplectic structure is still non-degenerate. We do this by introducing new phase space 
variables, which are similar to the usual non-relativistic centre of mass coordinates and momenta 
of a two particle system. The general solution to the equation P = M~y can obviously be written 

as 

p 1 = M 1 ~f -K~f(P), p 2 = M 27o + K~f((3). (1.8) 

The new variables are the energies of the particles, and K and (3 are polar coordinates defin- 
ing the spatial momentum and its direction. The rotating unit vector is introduced in ( |A.6| ). 
It defines the angular direction (3 in Minkowski space, thus in this case the direction of motion of 
the particles in the centre of mass frame. 

To express the mass shell constraints in terms of these variables, it is convenient to replace the 
mass parameters m x and m 2 by a total mass /i and a relative mass v, 

ii = m2 + mi, v = 77i2 — mi, < V < 11. (1-9) 

Without loss of generality, we assume that m 2 > mi. The special cases are v = 0, where both 
particles have the same mass, and v = \i, where at least one particle is massless. The same 
redefinition can be applied to the energy variables, 

M = M 2 + M 1 , V .U,-.U ; . \V\ < M, (1.10) 

where the inequality represents the positive energy condition for both particles. The mass shell 
constraints ( |1.2D are then given by 



1 2 8 2 2 8 

They can be simplified by taking the following linear combinations, 

V = C 2 -Ci = l —— 2 , £ = C 2 + d = K 2 + . (1.12) 

Let us solve the constraint V = for V, so that V becomes a simple function of M, 

What remains is a single mass shell constraint, which can be written as 

£ = K 2 -F(M)^0, where F(M) = (M ' ~ - ~ ^ . (1.14) 

We'll see later one that this is just a somewhat unusual way to write the familiar relation between 
the momentum and the energy of two relativistic point particles. The positive energy condition 
becomes a non-trivial condition to be imposed on M, namely 

\V\<M ^7<M & M>^. (1.15) 

ii M 
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Note that this implies M ^ 0, so that the various places where M appears in the denominator 
are not problematic. All together, the momentum vectors are parameterized by three independent 
variables M, K, (3. They have an immediate physical interpretation as the total energy, the spatial 
momentum, and the direction of motion of the particles in the centre of mass frame. 

The positions x k are subject to the constraints J = Sj and x\ = x\. The general solution 
can be parameterized by three other variables T, X, and 0, so that 

T = x° 1 = x° 2 , z = x 2 -x l = X 1 (<p). (1.16) 

The radial coordinate X > represents the relative position of the particles in space, and </> 
is the spatial orientation. It is also useful to introduce the relative position vector z, which is 
going to have a generalization for the interacting particles later on. The clock T represents the 
absolute time in the centre of mass frame. It is the 7 -coordinate of the equal time plane in 
figure 0. So far, it is an arbitrary function of the ADM time t, because we are still free to choose 
the parameterization of the world lines, and thus the labeling of the equal time planes by the 
unphysical coordinate t. 

Using all this, we can solve the equation J = S~( Q , and express the result in terms of the new 
phase space variables, 

M + V M - V 

xi = Tjo ~ x 2 = T 7o + _^_X 7 (0). (1.17) 

The relativistic version of the centre of mass frame is thus actually a centre of energy frame. At 
each moment of time, the centre of energy is located on the 7 -axis in figure 0. But nevertheless, 
let us stick to the notion centre of mass frame. The crucial point is that the positions of the 
particles are, like the momenta, specified by three independent variables. We have the clock T, 
the relative position X, and the orientation 0. 

To see that the six eliminated constraints were second class constraints, we have to compute 
the reduced symplectic potential. Inserting fll.8p and fll.l7p into fll.ip gives 



6 = K cos(/3 - 4>) dX + X K sin(/3 - 0) dcp - MdT. (1.18) 
To simplify this, we replace the polar momentum coordinates K and (3 by Cartesian coordinates 

Q = XK cos(/3-0), S = X K sin(/3 -(f)). (1.19) 
They define the radial and angular momentum. One can easily verify that S indeed satisfies 



S70, so it coincides with the previous definition. Moreover, inserting this into ( |1.8[ ) gives 



M — V g 7 (0) + gy(0) M + V Q 7 (0) + g 7 '(0) 
Pi = — 2 — 70 x ' P2 = — 2 — 70 X ' 

This tells us that Q/X is the component of the momentum parallel to the relative position, and 
S/X is the component orthogonal to it. This is the usual definition of a radial and angular 
momentum. If we use Q and S as the basic phase space variables, then the expression to be 
inserted into the mass shell constraint ( |1.14[ ) is 

= (1.21) 
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And finally, the symplectic potential simplifies to 

= X^QdX + Sd0- MdT. (1.22) 

This defines a non-degenerate symplectic structure Q — d©, and we read off the following 
non- vanishing Poisson brackets, 

{M,T} = -1, {Q,X} = X, {SA} = 1. (1.23) 

From this it is immediately obvious that M is the charge associated with time translations T i— > 
T — AT, and S is the charge associated with spatial rotations i— > + A<f>. 

What remains from the Hamiltonian is a single mass shell constraint £, and a multiplier £, 
which is some not further interesting linear combination of the original multipliers 

H = (£. (1.24) 

It is still an unphysical Hamiltonian. Its value is zero for physical states, and it generates the time 
evolution with respect to the unphysical ADM time t, which is formally a gauge transformation. 
Except for the reduced symmetry group, this looks very much like the Hamiltonian description 
of a single relativistic point particle. We have a six dimensional extended, or kinematical phase 
space 

V = { {M,Q,S;T,X,<f>) | X > 0, = + 2vr }, (1.25) 

and a single mass shell constraint 8, which defines the physical subspace and an associated gauge 
symmetry. It is this feature that we are going to exploit in the following, treating the two particle 
system in the centre of mass frame as if it was a single particle system. Of course, this is a well 
know way to solve a two particle system in non-relativistic mechanics, and we'll see that it also 
works for this almost trivial relativistic system. 

Complete reduction 

We can also go over to an effectively non-relativistic formulation, where no constraint and no 
gauge symmetry is left. We just have to impose another gauge condition. The most natural one 
is to require the ADM time t to coincide with the absolute time T in the centre of mass frame. 
This provides another constraint, and together with the mass shell constraint we get a new pair 
of second class constraints, 

K 2 = F(M), T = t. (1.26) 

Let us assume that the function F can be inverted. We can then define a four dimensional reduced 
phase space 

V = { (Q,S;X,<f>) | X>0, = + 2tt }. (1.27) 

It is just the usual non-relativistic phase space of the relative motion of two particles in a plane. 
To derive the symplectic structure and the Hamiltonian on?, we have to take into account that 
the constraints ( |1.26[ ) are explicitly time dependent. This implies a mixing of the Hamiltonian 
and the symplectic structure when we perform the reduction. We have to consider the extended 
symplectic potential on V, 

6-Hdt = X^Q dX + S dcj) - M dT - (£ dt. (1 .28) 
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To obtain the reduced structures on V, we insert the solutions to the equations K 2 = F(M) 
and T = t, and write the result as a combination of the reduced symplectic potential and the 
reduced Hamiltonian H, 

G-Hdt = X~ X Q dX + S d<f) - F-\K 2 ) dt. (1.29) 



Note that the last term in ( |1.28| ) vanishes, because the constraint £ = is now identically satis- 
fied. What comes out is 

= X- l QdX + Sd(f) {Q,X} = X, {S,(/>} = 1. (1.30) 

The Poisson brackets of the remaining variables are unchanged, and they are the usual non- 
relativistic ones. 

So, we now have an unconstrained Hamiltonian formulation of the free particle system, where 
the Hamiltonian represents the physical energy of the system, and generates the time evolution 
with respect to the absolute time in the centre of mass frame. Explicitly, one finds that 



H = F-\K 2 ) = ^K 2 + m x 2 + ^K 2 + m 2 2 . (1.31) 

Not surprisingly, this is just the relativistic energy of two point particles with the same spatial 
momentum K. The signs of the two square roots are fixed by the positive energy condition. A 
minus sign for one of the square roots would violate ( |1.15[ ), as one of the particles would then be 
in an antiparticle state. To express H as a function of the reduced phase space variables, we have 
to insert flOip . 

It is also quite instructive to look at the range of H . The right hand side of ( |1.31| ) is obviously 
minimal for K = 0, and it increases unboundedly with K. Thus, we have H > M min , where 

M min = mi + m 2 . (1-32) 

For K = the particles are at rest with respect to each other, and thus also with respect to 
the centre of mass. There is one exception, however, where such a state cannot be realized. If 



one of the particles is massless, then we have v — //, and the positive energy condition ( |1.15| ) 
requires that M > M min . In this case, the states with K = are excluded, and we have the 
stronger condition H > M min . Clearly, this is because a massless particles cannot be at rest, and 
consequently a state with vanishing momentum K does not exist. 

Let us also consider the non-relativistic limit. For small momenta K, the energy H either starts 
off linearly or quadratically with K, depending on whether a massless particle is present or not. 
Let us consider the case where both particles are massive. If we then expand the right hand side 



of fll.31p up to the second order in K, we get 



H w M min + — K 2 , where m = mim2 . (1.33) 
2m mi + W2 

This is the usual non-relativistic relation between the spatial momentum K and the energy H in 
the centre of mass frame. The parameter m is the reduced mass of the two particles system. The 
difference between the total energy H and the total rest mass M min is the non-relativistic kinetic 
energy. For large momenta K, on the other hand, we find the usual relativistic behaviour H w 
2K. The energy becomes a linear function of the spatial momentum. The factor of two arises 
because we have two particles with the same momentum. For some typical mass parameters, the 
relations between K and H are shown as the broken lines in figure |[ 
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Trajectories 



There are now two alternative ways to describe the two particle system. We can either use the 
constrained Hamiltonian formulation based on the extended phase space V . Or we may use the 
unconstrained formulation, based on the reduced phase space V. Consider first the constrained 
formulation. In this case we have six independent phase space variables, and the Hamiltonian is 
given by H = (£, where £ is the mass shell constraint and ( is an arbitrarily chosen function of 
the unphysical ADM time t. It is not difficult to derive the resulting time evolution equations and 
to solve them. 

The energy M and the angular momentum S are of course preserved charges. The same holds 
for the spatial momentum K, as it is only this combination of the phase space variables Q, S, 
and X that enters the Hamiltonian. Thus, 

M = {H,M} = 0, S = {H,S} = 0. K = {H, K} = 0. (1.34) 

For the relative position X and the conjugate radial momentum Q, we find 

X = {H,X} = 2Cj, Q = {H,Q} = 2C^^ = 2(K 2 . (1.35) 
And finally, the brackets of H with and T are given by 



{H,4>} = 2(j^, T = {H,T} = (F'{M). (1.36) 



It is clear that the same equations of motion arise in the unconstraint formulation, except that 
M and T are not independent variables, H is replaced by H, and instead of the multiplier £ 
the function 1/F'(H) appears, which formally implies that T — 1, which is consistent with the 
gauge condition T = t. 

We can easily solve these differential equations step by step. Those for M and S are trivial, 
stating that M(t) = M and S(t) = S for some constants M and S . Moreover, K(t) = K 
is also a constant of motion. It is related to M by the constraint Kq = F(M ). Thus K and 
M are not independent, and they are subject to the previously derived restrictions. For massive 
particles, we have M > M min and K > 0. If at least one massless particle is present, then we 
have the stronger condition M > M min and K > 0. 

Using this, it is easy to solve the equation of motion for Q. The general solution is Q(t) = 
Qo + e(t)KQ, where Qo is some integration constant. The function e(t) is determined up to a 
constant by e(t) = 2((t). For K > we can obviously choose this free constant so that Q = 0. 
On the other hand, K = implies, by definition fll.21| ), that both Q and S must be zero. In this 
case we also have Q = 0. Therefore, the general solutions found so far are 

M(t) = M , S(t) = S , K(t)=K , Q{t) = e{t)Kl (1.37) 

If we insert this into the time evolution equations for X and 0, then we can integrate them. The 
result is 

X(t) = ^Rl + e 2 (t)Kl (f>(t) = 0o + arctan (e(t) , (1.38) 
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where R > and </> are two more integration constants. And finally, we can also solve the 
equation of motion for T, which gives 

T(t)=T + e(t)^^. (1.39) 

These are the most general solutions to the time evolution equations on V, provided by the 
Hamiltonian H = (£. They are parameterized by an arbitrary function e(t), representing the 
gauge freedom, and six integration constants M , K , S , T , Ro, (j>o- Only four of them are 
independent. The constraint relates the momentum K to the energy M , and additionally there 
is a relation between S , X and K , which follows from the definition Q1.21D , 

K 2 =F(M ), K R = \S \. (1.40) 

There are two ways to look at these solutions. If we keep t fixed and vary e(t), then we pass 
along a gauge orbit generated by the constraint £. If we instead fix the function e(t) and vary t, 
then we see a trajectory in the phase space V, which describes the time evolution generated by 
the Hamiltonian H in a particular gauge. 

The derivation of the trajectories on the reduced phase space V yields the same result, with the 
same integration constants, however with the gauge function replaced by 

e(t) = 2 - ~ T ° , (1.41) 

w F'{M ) 

which obviously implies T(t) = t. The inverse of F'(M ) is always well defined because, within 
the range of M allowed by the positive energy condition, the function F(M) is monotonically 
increasing with M. So, the trajectories on V are parameterized by the same four independent 
integration constants, but we no longer have any gauge freedom. 

Let us now look a little bit closer at the various trajectories. Consider first the case where 
So ^ 0. The second relation in fll.40p then implies that both K > and R > 0, and Q1.38D 



describes a straight line in polar coordinates. The particles approach each other from a direction 
0i n . At the time T they reach a minimal distance R with spatial orientation (fi . And finally 
they separate again, moving into a spatial direction out . The directions in and <p ont are always 
antipodal, and both are orthogonal to 4>$, 

IT 7T 

4>m = 0o-2 sgn So, 4>out = 0o+2 sgn S <j) ut = 4>m + 7T sgn S Q . (1.42) 

The sign of the angular momentum S tells us on which side the particles pass each other. This is a 
trivial scattering process, because the directions in and out always differ by tt. Nevertheless, the 
scattering picture is quite useful. There will by a real scattering when the gravitational interaction 
is switched on. 

The constant of motion K = K has a useful interpretation in this picture. It represents 
the momentum which is canonically conjugate to the distance between the particles, when this 
distance is very large. More precisely, consider the pair of canonically conjugate variables 

P = Q/X, R = X {P,R} = 1, (1.43) 
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replacing Q and X. Thus, R also represents the distance between the particles in space, but P 
is now the momentum which is canonically conjugate to it, in contrast to Q which is not exactly 
conjugate. For the trajectories derived above, we have 

R(t) = ^Rt + e*(t)Kl P(t) = £(t) f (1.44) 

^Rl + e\t)Kl 

In the limit t — > ±00, or equivalently e(t) — > ±00, we have R(t) — > 00 and P(t) — > ±K . The 
constant of motion K = K is equal to the momentum P when the particles are far apart. This 
will be useful to know when the system is quantized later on. The eigenvalue of K is then the 
inverse radial wavelength of the wave function at spatial infinity. At the classical level, K is the 
momentum of the incoming and outgoing particles in the scattering process. 

A special situation arises when Sq = 0. Then, according to fll.4(jp , we have either K = or 



Ro = 0. The first case, Sq = 0, K = 0, and R > 0, is the trivial one. In Q1.38D , we may replace 
K%/ 'So by Sq/ Rq. Then we get the simple trajectory X(t) = R and 4>(t) = O - The particles 
are at rest, at a spatial distance Rq and with angular orientation O - The integration constant T 
is in this case redundant, and the energy M = M min is equal to the sum of the rest masses. Of 
course, such a trajectory only exists when both particles are massive, as otherwise we must have 
K Q > 0, because massless particles cannot be at rest. 

The more interesting case is S = 0, R = 0, and K > 0. The argument of the arctan is then 
ill defined. But we can still consider the trajectory as a limit. In the limit S — ► and R — > 0, 
with Kq and all other integration constants fixed, we get 

X{t) = \e{t)\K , (Pit) = O + I sgnSo sgne(f). (1.45) 

The crucial point is now that the limit still depends on sgn S , thus on the direction of the limit. 
But if we change the definition of the integration constants slightly, we may also write 

«*>={£. tltt «- 46 > 

where the sign no longer matters. Clearly, this is just the definition of a straight line through 
the origin of the polar coordinate system. The particles approach each other from a direction in , 
touch each other at t = T without interaction, and separate again into the opposite direction out . 
The case Sq = is only special because of a coordinate singularity of the phase space variables. 
We could avoid this by introducing a globally well defined chart on V and V, replacing the polar 
coordinates X and by Cartesian coordinates, and the radial and angular momenta Q and S by 
the appropriate conjugate variables. 



Schrodinger quantization 

Let us now, finally, perform the quantization of the free particle system. The most straightforward 
way is to start from the reduced classical phase space V = {(Q, S; X, 0)} with symplectic 
potential Q1.30D , 

6 = X^QdX + Sd(j). (1.47) 
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This is an unconstrained formulation. We can apply the standard Schrodinger quantization pro- 
cedure. We use a position representation, where the wave function is given by ip(x,(p), with x 
and cp being the eigenvalues of the operators X and 0, 

Xip(x,(p) = xif)(x,(p), <j)t/)(x,<p) = ipip(x,ip), (1.48) 

with x > and ip = (p + 2n. Actually, the second equation is of course only well defined with 
replaced by a periodic function of 0, as only those functions are well defined on the classical 
phase space. Concerning the coordinate singularity at x = 0, we shall follow the concept in 
It can later also be applied to the interacting system without any essential modification. We 
impose a periodicity condition on the wave function, 

ip(x,cp + 2ir) = e 27riA ^(x,^), (1.49) 

and we require the wave function to be finite at x = 0. The parameter A is some fixed real 
number. It represents a quantization ambiguity that typically arises when the classical phase 



space is not simply connected [p3||. A special case occurs when the two particles have the same 
mass. They can then be regarded as identical, and replacing ip with ip + ir already takes us back 
to the same state. In this case, we have the stronger relation 

ip(x,ip + ir) =e 7TiX i;(x,ip), (1.50) 

and the parameter A defines the statistics of the particles. If A is an even integer, then the particles 
are called bosons, and for odd integers they are called fermions. Note that this only refers to the 
statistics, not to the internal structure of the particles. There is no spin statistics theorem for this 
simple toy model. 

For non-integer values of A, the particles are usually called anyons [fTH|. Anyons can only arise 
in three spacetime dimensions. In higher dimensions, the fundamental group of the phase space 
of two identical particles is Z 2 rather than Z, which implies that A has to be an integer, which 
is even for bosons and odd for fermions. Another special feature in three spacetime dimensions 
is that the non-trivial phase factor in (|1 .49[) also shows up when the particles are not identical. 



Slightly abusing the language, we shall then also refer to A as the statistics parameter. We should 
also note that for identical particles the real number A is defined modulo two, otherwise modulo 
one. 

For ip(x, (p) to represent the usual probability amplitude in polar coordinates, we define the 
scalar product to be 

(^1^2) = J xdxd(pi> 1 (x,(p)i/j2(x,(p). (1-51) 
The momentum operators are defined so that the commutators are — ih times the Poisson brackets, 

[Q,X]=-ihX, [§,$] = -ih. (1.52) 

They become self-adjoint operators if we set 

d ~ d 

Qip(x,ip) = -ih — xip(x,<p), Sip(x,<p) = -ih — ip{x,cp). (1.53) 

OX Oip 
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Note the extra factor of x in the operator representation of Q. It has to appear in the given or- 
dering, since otherwise the resulting operator is not self-adjoint with respect to the inner product 
( |1.5ip . We should also note that the operators for Q and S are uniquely defined by the commu- 
tation relation, up to redefinitions that can be compensated either by a phase transformation of 
the wave function, or by redefining the periodicity condition ( |1.49[ ). In other words, the statistics 
parameter A is the only relevant quantization ambiguity. 

To derive the energy eigenstates let us first consider the operator K 2 , of which the Hamiltonian 
H = F~ l (K 2 ) is a function. We choose the Hermitian ordering 

K 2 = X^iQ 2 + S 2 )X-\ (1.54) 

In the representation ( |1.53[ ), K 2 is then proportional to the usual Laplacian in the polar coordi- 
nates, 

K 2 ib(x,tp) = —h 2 Aib(x,(p), where A = x~ 1 — x— — \-x~ 2 —- — . (1.55) 

ox ox dip dip 

To make K 2 self-adjoint, a boundary condition at x = is needed. Following $23j\, we choose 
this condition so that the wave function remains finite at x = 0. The normalized eigenfunctions 
are then parametrised by two quantum numbers, k and s, and given by the Bessel functions of 
the first kind, 

x (k, 8 ;x,<p) = ^=e im fJ\ a \(kx). (1.56) 
y An 

They are normalized so that 

J xdxdipxiki, sr,x,ip) x(k 2 , s 2 ',x,<p) = k~ l 5(k 2 - h) 5 S2 - Sl . (1-57) 

As expected, the quantum number k is the inverse radial wavelength of the wave function at 
infinity. We already expected the eigenvalues of K to have this property, because for large 
distances of the particles the momentum K is equal to the canonical radial momentum P = 
Q/X. The quantum number s is the inverse angular wavelength, and of course it represents 
the eigenvalue of the angular momentum. The possible values of s can be inferred from the 



periodicity condition ( |1.49[ ). It implies 

e is( V +2n) = e 2niX ^ ^ S G A + Z. (1.58) 



For identical particles, the stronger condition ( |1.50D implies 



e «(¥>+T) = e «iA e w ^ s G A + 2Z. (1.59) 

It follows that the eigenvalues of K and S are 

K — hk, S = hs. (1.60) 

The spectrum of K is continuous and positive, and S is quantized in steps of h, or 2% for identical 
particles. The actual values of S are determined by A. In the special case of two identical 
bosons, for example, the total angular momentum is an even multiple of h, and for two identical 
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fermions it is an odd multiple of h. This is what we usually find for bosons and fermions in 
higher dimensions as well. The peculiar feature of anyons is that the spectrum of the angular 
momentum is shifted by a non-integer multiple Xh. 
The completeness relation of the eigenfunctions reads 

/ kdkx{k,s]Xi,tpi) x(k } s;x 2 ,(f2) = x' 1 5(x 2 - x x ) 5\(ip 2 - <pi), (1-61) 

s •* 

where s takes the values given above, and 5\ is a periodic delta function. It fulfills the same 
periodicity condition as the wave functions, hence ( |1.49| ) or Q1.50D , depending on whether the 
particles are identical or not. 

The momentum eigenstates x{k, s; x, (p) are also the energy eigenstates. The Hamiltonian is a 
function H = F~ l (K 2 ) of K, so that the energy becomes a function of the quantum number k. 
It is useful to introduce a dispersion relation, which represents the classical relation between the 
momentum K and the energy H as a quantum relation between the inverse radial wavelength k 
at infinity and the frequency to of the wave function, 

u(k) = h^ 1 F^ 1 (h 2 k 2 ). H X {k,s;x,ip) =hu{k)x{k,s;x,ip). (1.62) 

It follows that the spectrum of H coincides with its classical range Af min < H < oo. And finally, 
we can write down the general solution to the time dependent Schrodinger equation, which is a 
superposition of energy eigenstates with appropriate frequencies, 

ip(t;x,ip) = J2 I kdke~ i ^ k)t i/;(k,s)x(k,s;x,if). (1.63) 

The time-independent function ip(k, s) represents the probability amplitude in momentum space. 
Clearly, this is a general superposition of plane waves, written down in polar coordinates using 
the Bessel functions. The Schrodinger quantization of the free particles system is completely 
straightforward and gives the expected result. 

Dirac quantization 

There is not much more to be said about this very simple dynamical system. However, in case 
of the Kepler system it turns out that a simple Schrodinger quantization like this is not possi- 
ble. The reasons are not immediately obvious at this point, so let us not discuss them here. 
Instead, let us look for an alternative method that leads to the same result. A possible alter- 
native quantization is the Dirac procedure, based on the six dimensional extended phase space 
V = {(M, Q, S; T, X, 4>)}, with symplectic potential flL22D , thus 



= X-'QdX + Sd4>- MdT. (1.64) 

To quantize this phase space, we choose a wave function tp(r, x, cp), which additionally depends 
on the eigenvalue r of the clock T. In addition to the operators ( |1.48D and ( |1.53D we have 

d 

Tip(r,x,tp) = Til)(T,x,ip), Mip(r,x,(p) = ih— ip(r,x,(p). (1.65) 

or 
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All the basic operators are then self adjoint with respect to the scalar product 

(i)i\i>2) = / dr xdxdpipi(r, x, p)ip 2 (r, x, <p). (1.66) 



On this extended Hilbert space, we have to impose the constraint £ = K 2 — F(M), and the 
positive energy condition M > sfjw. The constraint is given as a function of K and M, so it is 
useful to diagonalize these operators first. Clearly, the eigenstates are again given by the same 
Bessel function, 

X{uj,k,s;r,x,cp) = — j= e~ iuJT e istp J\ s \{kx). (1.67) 



/2tt 

The normalization now reads 



dr x dx dip fci, Si; r, x, <p) x(^ 2 , h, s 2 ; r, x, p) = 

2ti5{u 2 - w x ) k- 1 5{k 2 - fci) 6 aa - ai . (1.68) 

The spectrum of K and S is the same as before, and that of M is real and continuous. The 
eigenvalues are 

M = hu, K = hk, S = hs. (1.69) 

It is then very easy to solve the constraint. The states which are annihilated by the operator 
£ = K 2 — F(M) and satisfy the positive energy condition M > ^fjw are those where the 
quantum numbers are related by the dispersion relation cu = u(k). A general physical state is 
given by 

ip{r,x,(p) =53/ kdkip{k,s)x(u(k),k,s;T,x,<p), (1.70) 

s J 

where ip(k, s) is again the wave function in momentum space. It is useful to split off a radial 
wave function and write the result as 

V>(t,x,</?) = 53 f kdke- [ ^ T e lsv ^{k,s)C{k,s;x), (1.71) 

s J 

where 

C(k,s;x) = -^=J\.\(kx). (1.72) 

We shall later compare this radial wave function to that of the Kepler system, and from this we 
will be able to read off the basic effects of quantum gravity in this toy model. 

Formally, the physical wave functions ( |1 .70[ ) in the Dirac approach are exactly the same as 
Ql .63[ ) in the Schrodinger approach. There is only the following conceptual difference. The wave 



function in the Dirac approach is not a solution to the time dependent Schrodinger equation, but 
a solution to a generalized Klein Gordon equation, thus a constraint equation. In the Schrodinger 
formulation, the wave function depends by definition on the physical time t, whereas in the Dirac 
formulation, the time dependence is encoded implicitly in the dependence of the wave function 
on the eigenvalue r of the clock, which is one of the classical phase space variables. 

In the Dirac approach, there is no time evolution of the state with respect to the ADM time t, 
because this is an unphysical coordinate. The Hamiltonian annihilates all physical states, and the 
Schrodinger equation is void. To say something about the physical time evolution in the Dirac 
quantization, we have to refer to an operator that represents a clock. The free particle system is 
in this context a nice toy model to explain how the problem of time arises in quantum gravity, 



and how it can be solved by considering operators representing physical clocks [24|. 
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2 The classical Kepler system 



Now we are going to switch on the gravitational interaction. As already mentioned in the in- 
troduction, we shall first consider an alternative way to construct the Kepler spacetime, which 
includes the definition of a centre of mass frame, based on the asymptotic structure of the space- 
time far away from the particles. We shall thereby focus on the geometrical aspects, and skip all 
technical details and proofs. They can be found in pCp. We shall then introduce a set of phase 
space variables. They provide position and momentum coordinates, similar to those of the free 
particles. But at the same time they also specify the geometry of space at a moment of time, thus 
providing the of ADM variables of general relativity. 

On the phase space spanned by these variables, we can then introduce a Hamiltonian and a 
symplectic potential. The former can be inferred from some straightforward geometric consid- 
erations. The latter can only be derived from the full theory of Einstein gravity as a field theory. 
As this is not part of this article, we can here only refer to the general derivation of the phase 
space structures for a multi particle model in [J]~0|] . On the given phase space, we can then derive 
and solve the equations of motion, and finally we shall briefly discuss the various trajectories and 
compare them to the free particle trajectories. 

The static cone 

To explain the definition of the centre of mass frame, and to show how the phase space variables 
of the Kepler system specify the geometry of space, it is useful first to consider a simplified 
example. Let us assume that the particles are at rest with respect to each other, so that the 
spacetime is static. It is then possible to introduce a globally defined absolute time coordinate 
T, and the spacetime can be foliated by the surfaces of constant T = T. In the following, a 
letter with a tilde always denotes a coordinate on the spacetime, whereas all other symbols are 
configuration or phase space variables, or functions thereof. 

In the static spacetime, the surface T = T represents an instant of time, like the equal time 
planes for the free particle system in figure [Z[ All slices have the same geometry, and the space- 
time is just the direct product of a fixed space with a real line. The Einstein equations require this 
space to be locally flat, and there must be two conical singularities, representing the particles. In 
other words, the space is a conical surface with two tips. We denote them by n k , with k — 1, 2. 
The deficit angles at the tips are equal to 8%Gmk. What we would like to find is a suitable set of 
configuration variables, describing the geometry of such a surface. And we are also looking for 
an appropriate definition of a centre of mass frame. 

For this purpose, we embed the conical surface in a special way into a Euclidean plane, as 
shown in figure 0. The first step is to introduce a geodesic A, which connects the two tips 7Ti 
and 7r 2 on the conical surface. There is always a unique such geodesic, because apart from the 
two tips there is no curvature. When the conical space is cut along the line A, it becomes a 
surface with a hole in the middle. The points nk are located on a circular boundary. As shown in 
figure ^|(a), the resulting surface is a subset of an ordinary cone, with a certain region around the 
tip taken away. The total deficit angle of the big cone is the sum of the deficit angles at the tips 

7T fc , thus 8^(7(772! + m 2 ). 

Let us introduce on the big cone a polar coordinate system (R, 0), so that R defines the radial 
distance from the tip, and <p is an angular coordinate with a period of 2n. Including the absolute 
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Figure 3: A double cone with two tips can be constructed by cutting off the tip from a single 
cone. The cut lines A, which are afterwards identified, are two geodesies connecting the 
points 7Ti and 7r 2 . If the spacetime is static, then it is the direct product of this space with 
a real line. If the conical surface is cut not only along the geodesic A, but also along two 
radial lines 771 and r/ 2 , then the two half spaces A± can be embedded symmetrically into a 
Euclidean plane. 



time T as a third coordinate, the spacetime metric can be written as 

ds 2 = -df 2 + dR 2 + (1 - AG{m x + m 2 )) 2 R 2 d0 2 . (2.1) 

This defines a static cone with a total deficit angle of 87rG(mi + m 2 ). The conical surface in 
figure |3|(a) becomes a surface of constant T = T in this spacetime. Roughly speaking, we can 
say that the static Kepler spacetime is a cone with a tip cut off, and with the cut lines identified. 

From the physical point of view, the situation is as follows. Far away from the particles, say 
at spatial infinity, the spacetime looks like the gravitational field of a single particle with mass 
mi + m 2 . It is reasonable to identify the rest frame of this fictitious particle with the centre of 
mass frame of the universe. The centre of mass frame is thus defined by a certain property of the 
spacetime metric at spatial infinity, which does not directly refer to the particles. It is a kind of 
asymptotical flatness, which is explicitly defined in terms of the conical coordinates (T, R, 0). 

The fictitious tip of the big cone in figure |3|(a) is the apparent location of the centre of mass, as 
seen by an observer at infinity. In analogy to the free particle system, let us identify the centre of 
mass frame with the reference frame of such an external observer. With respect to this reference 
frame, we define the absolute positions of the particles in spacetime as the conical coordinates 
(T k , R k , <p k ) of the points ix k on the big cone. Again in analogy to the free particles, these absolute 
coordinates are not independent. They are specified by three independent relative coordinates. 

At each moment of time, both points n k are located on the same surface of constant T = T. 
Moreover, for the cut line A to be of the same length on both sides of the big cone, the points ir k 
must be located on two opposite, or antipodal radial lines. Hence, we must have 

T 1= T, T 2 = T, 1 = 0±7r, 2 = 0. (2.2) 
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The variables T and <p are straightforward generalizations of the corresponding free particle vari- 
ables. The former defines a clock, the latter represents the angular orientation of the particles 
with respect to the reference frame. As a third independent coordinate, we introduce the distance 
R between the particles. It is the length of the geodesic A. It is a simple exercise in conical geom- 
etry to compute this. The result can be expressed as a function of the conical radial coordinates 
R k and the total deficit angle, 

R 2 = R\ + R\ + 2 R x R 2 cos(47rG(ma + m 2 )). (2.3) 

Finally, there is one more consistency condition. The deficit angles at the points ir k must be equal 
to 8nGm k . So far we have only fixed the total deficit angle of the big cone, which is the sum of 
the two individual deficit angles. We still have to fix the way this deficit angle is distributed over 
the two tips. This yields another relation between the radial coordinates R k . It is again a simple 
exercise in conical geometry to show that 

sin(47rGm 2 ) sin(47rGmi) 

Ri = — — 7- — —, rr R R2 — ~. — ~r. — 7-77 77 R- (2.4) 

sm(47rG(mi + 1112)) sm(47rG(mi + m2)) 

So, what is the conclusion? There are obvious three independent variables T, R, and (p. They 
define the geometry of the space manifold at a moment of time, and the way this space manifold 
is embedded into the spacetime manifold. In this sense, they provide some kind of discretized 
ADM variables of general relativity. On the other hand, from the particle point of view, the 
same variables define the relative position of the particles with respect to each other, as well as 
their absolute positions with respect to the centre of mass frame. The definition of the relative 
coordinates is very similar to (|1.16[) for the free particles. 



We cannot say anything about the momentum variables at this point, because we are only 
considering the static states. We still have to generalize these concepts for moving particle. But 
let us stick to the static case for a moment, and let us introduce a slightly different representation 
of the same conical geometry of space. Let us not only cut the space manifold along the geodesic 
A, but also along two antipodal radial lines r] k , extending from the points 11 k to infinity. Hence, i] k 
is a line of constant T = T k and <p = 4>k in the big cone, which lies inside the surface of constant 
conical time T — T. The conical surface is then divided symmetrically into two half spaces, 
which we denote by A + and A_. 

Each half space A± is bounded by three edges, denoted by r)x±, A±, and r]2±, and it has two 
corners, called m± and 7r 2 ±. The two half spaces are flat and simply connected. They can be 
embedded into a Euclidean plane in the following unique way, which is shown in figure |3|(b). 
The apparent position of the fictitious centre of mass, hence the tip of the big cone, is mapped 
onto the origin of the plane. Moreover, the edges A± are mapped onto two parallel straight lines, 
with angular direction <p. This fixes the embedding of the half spaces into the plane completely. 

To describe the embedding more explicitly, it is convenient to think of the Euclidean plane as 
an equal time plane in Minkowski space. We choose it to be the plane with 70-coordinate T. The 
geometry of the two half spaces can then be specified as follows. The edges A± are represented 
by two spacelike Minkowski vectors z±. They are given by 

z ± = R 1 {cf ) ). (2.5) 
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The end points of the edges A± are the corners Hk±- They are located at the following points in 
Minkowski space, 

x 1± = T-y -R 1 7(0 T 47rGmi), x 2± = T lo + R 2l (<P ± 4nGm 2 ), (2.6) 

where the coordinates Rk given by (|2~4l) . One can easily verify that this implies z± = x 2 ± — xi±. 
And finally, the edges 7]k± are radial lines, extending from the corners 7Tfc± to infinity. A radial 
line in Minkowski space is a geodesic which is orthogonal to the 70-axis, and intersects with it 
if it is extended beyond its end point. 

The advantage of figure |3](b) is that it uses less dimensions to visualize the conical geometry 
of space. In contrast to figure ^](a), we do not need an auxiliary third Euclidean dimension to 
embed the big cone. This will be useful when we now consider the moving particles, because 
then the space is no longer flat. We can also read off the deficit angles of the particles directly 
from figure |3|(b). The deficit angle of the particle iti is the angle between the edges r)i+ and 771- 
in the plane. According to (^!6|), this angle is equal to SitGmi. The deficit angle of the particle 
tt 2 is the angle between the edges 7/ 2 _ and r] 2+ on the other side, which is equal to 8irGm 2 . 

The disadvantage of figure |](b) is that the definition of the centre of mass frame is less obvious. 
However, it is actually only the relation between the geometry of space and the centre of mass 
frame, which we need in the following. And this can still be read off from figure [3](b). First, 
we need to know which time slice is represented by the given space. This is specified by the 
clock T, and according to ( [ZIjI ) this is also the 7 -coordinate of the plane shown figure |3](b). And 
secondly, we need to know the spatial orientation of the particles with respect to the reference 
frame. This is specified by the symmetry axis of figure |3|(b). It is the solid line in the middle, 
which is a radial line with angular direction 7(0). 



The spinning cone 

When the particles are moving, we have to generalize two aspects of this construction. The 
geometry of space is no longer constant in time. The previously introduced relative position co- 
ordinates become time dependent, and we also have to introduce conjugate momentum variables. 
The time evolution will finally be provided by some Hamiltonian, which we still have to derive. 



Moreover, we also have to replace the conical metric ([Zj]) of the big cone by a more general 
conical geometry. Let us first consider this second aspect. 

If the particles are in motion with respect to each other, then the fictitious centre of mass 
particle receives a variable mass M and a variable spin S, where M is the total energy and S is 
the total angular momentum of the system. The gravitational field of such a particle is a spinning 
cone [0]. Using the same conical coordinates as before, the metric of a spinning cone can be 
written as 

ds 2 = -(df + AGS d0) 2 + dR 2 + (1 - AGM) 2 R 2 d0 2 . (2.7) 

This is obviously a generalized version of Q2.1|). The previously considered static case is recov- 
ered if we set M = mi + m 2 and S — 0. But nevertheless, we can still use the conical coordinates 
(T, R, 4>) to define the centre of mass frame of the universe, and we may also identify this with 
the reference frame of some external observer sitting at infinity. 

To visualize the geometry of the spinning cone, let us consider a surface of constant conical 
time T = T. We can still say that each such surface represents an instant of time in the centre of 
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mass frame. Moreover, all such surfaces have the same geometry, as long as we do not insert the 
moving particles. The spinning cone is no longer static, but still stationary. A surface of constant 
conical time is not flat, but it can locally be embedded into Minkowski space. The following 
local isometry maps a segment of the spinning cone into Minkowski space. We fix some angular 
direction O and define 

(T,R,$) i- (T + 4G5(0-0 O )) 70 + i?7(0 - 4GM (0 - O )). (2.8) 

One can easily check that this is indeed an isometry. It is only a local isometry, because it does 
not respect the periodicity of 0. A surface of constant T = T of the spinning cone is mapped 
onto a screw surface in Minkowski space. 

A screw surface in Minkowski space is defined by a family radial lines, where the 7 - 
coordinate increases linearly with the angular direction. The slope of this screw surface is deter- 
mined by the parameters M and 5. If we pass once around the spinning cone, thus if we increase 
by 27r, then the total angle covered by the screw surface in Minkowski space is 2n — 8ttGM, 
and the amount that we move forward in time is 8nGS. The spinning cone has a deficit angle of 
8nGM and a time offset of 8irGS. In units where the velocity of light is one, the latter is a length 
or time if 5 has the dimension of an angular momentum. 

The central axis of the spinning cone at R = is mapped onto the 7 -axis in Minkowski space. 
This is the apparent world line of the centre of mass, as seen by an observer at infinity. There is 
a critical radius R = AGS/ (1 — AGM). For R < R , the screw surface is timelike, at R — Rq it 
is lightlike, and for R > Ro it is spacelike. There are thus closed timelike curves in the spinning 
cone, passing through the region R < R . When we insert the particles, this critical region 
will be cut away, so that the actual Kepler spacetime has no closed timelike curves. Outside the 
critical region the causal structure of the spinning cone is well defined. 

To insert the particles, we have to generalize the construction in figure ^fb). Remember that 
the plane in which this construction took place was an equal time plane in Minkowski space, with 
7o-coordinate T. And remember also that the symmetry axis was a radial line in this plane, with 
angular direction 7(0). We shall now lift this picture into the third dimension. The half spaces 
A± are tilted in a certain way, while the symmetry of the figure is preserved. The deformation is 
so that finally the half spaces become two segments of the previously considered screw surface. 
Hence, they fit into the spinning cone as a surface of constant conical time T — T, 

As indicated in figure |], the deformation is as follows. The corners n k ± are raised and lowered 
symmetrically by 2nGS. Moreover, the mass parameters rrik appearing in (|2.6|), which define 



the deficit angles of the particles, are replaced by two independent energy variables Mfc. This 
is because for a moving particle the deficit angle is actually not proportional to its mass, but 
proportional to its energy [O]. We shall look at this more closely in a moment. As we want the 
symmetry axis of the figure to be preserved, the locations of the corners iT}.± are then given by 

x l± = (T ± 2ttGS) 70 - Ri 7(0 =F AttGMi), 

x 2± = (T t 2nGS) 7o + R 2 7(0 ± 4ttGM 2 ). (2.9) 



The radial coordinates R^ are still to be determined. We recover the static case ( |2.6[ ) if we set 
Mk = m k and 5 = 0. The edges rjk± are still defined as the unique radial lines in Minkowski 
space which pass through the corners iik±. As they are now located at different time levels in the 



24 



Figure 4: To define the geometry of space at a moment of time if the particles are not at 
rest, the construction in figure |3|(b) has to be lifted into the third dimension in the embedding 
Minkowski space. The half spaces are then two segments of a screw surface. This surface 
represents an instant of absolute time in the centre of mass frame. The geometry of the two 
half spaces, the way they are embedded into the reference frame, and also the locations of 
the particles in space are specified by a set of six independent phase space variables. 



embedding Minkowski space, each half space becomes a segment of a screw surface. The time 
offset of each half space is AtcGS, and the angle in Minkowski space covered by each half space 
is 7i - AnGM, where M = M l + M 2 . 

Let us glue the two half spaces together along the edges r] k ±. The result is a screw surface 
with a deficit angle of 8ttGM and a time offset of 8nGS. We can still use figure |](a) as a 
schematic representation of this surface. It is a conical surface with a hole in the middle or, 
roughly speaking, a cone with the tip cut off. It can be embedded into a spinning cone, as a 
surface of constant conical time T = T. Thus it represents an instant of time in the centre 
of mass frame, and the clock T represents the absolute time. The angular orientation of this 
embedding is fixed in the same way as before. The geodesies 77^ are two antipodal radial lines 
with conical coordinates T = T k and <p = 4>k, where 

T 1 = T, T 2 = T, 0i = 0±vr, 2 = 0. (2.10) 

The edges A± are the geodesies in Minkowski space, connecting the corners Hk±- As these are in 
general not contained in the screw surfaces, we actually have to deform the half spaces slightly 
in the neighbourhood of these edges. But this does not affect the embedding into the spinning 
cone. We can still assign two Minkowski vectors z± to the edges A±. They must be spacelike for 
the half spaces to be spacelike surfaces. And they can still be chosen so that both have the same 
direction <p, parallel to the symmetry axis. Moreover, we already know that the 70-components 
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of these vectors must be equal to the time offsets 47rG5' of the two half spaces. Hence, we have 



z± = ^JR 2 + (AirGSy 7 (0) =f AttGS 7o , (2.11) 

where R is again the distance between the particles. On the other hand, also have z± = x 2 ± — 
Xi±, which implies that the radial coordinates R k are given by the following generalization of 

(M), 



sin(47rG(Mi + M 2 )) v v ; 
sin(47rGM 1 



So, after this somewhat technical construction, what did we get? There are now six independent 
configuration variables Mi, M 2 , S, T, R, and 0. From the point of view of general relativity, 
they specify the geometry of space at a moment of time, providing a kind of discretized ADM 
variables. From the particle point of view, they define both the relative position of the particles 
with respect to each other, as well as their absolute positions with respect to the centre of mass 
frame. 

The relative position of the particles is specified by the geodesic distance R and the angular 
orientation 0. The absolute positions are given by the conical coordinates (T k , R k , fc ) of the 
points 7r fc in the spinning cone. All this is analogous to the previously considered static case, and 
also to the free particles, where the relative coordinates are X (or R) and 0, and the absolute 
positions are given by ( |1.17[ ). What makes the definitions a little bit more complicated here 



is that the reference frame is a conical frame, and not a Minkowski frame. What we do not 
see so far is in which sense the energy and momentum variables are conjugate to the position 
coordinates, how they evolve in time, and what is also still missing is a variable that replaces the 
radial momentum Q. 



Transition functions 

Let us have a closer look at the way the edges in figure [| are glued together. The embedding 
Minkowski space can locally be identified with the spacetime manifold, if the world lines are 
excluded. More precisely, it provides an atlas of the spacetime manifold in the neighbourhood 
of the embedded space manifold at the given moment of time. The atlas consists of two charts, 
and each chart contains one half space as a spacelike surface. There are three overlap regions 
between the two charts, corresponding to the edges r)i±, \±, and r/ 2 ±. The transition functions in 
these overlap regions are isometries of the embedding Minkowski space, hence Poincare trans- 
formations. 

Since in the following only the rotational components are of interest, let us ignore the trans- 
lations that are involved. We can then say that the edge ijx- is mapped onto the edge r/i + by a 
certain Lorentz rotation. It can be represented by an element hi E SL(2) of the spinor represen- 
tation of the three dimensional Lorentz group. Similarly, the edge r^ 2 is also mapped onto the 
edge 7/ 2 + by some Lorentz rotation, which is given by h 2 E SL(2). And finally, A_ is mapped 
onto A + by yet another Lorentz rotation g E SL(2). The transition functions tell us how the half 
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spaces are to be glued together along the edges. The question is thus, are they already determined 
by the variables introduced so far, or is there still some ambiguity? 

Consider first the edges r/i_ and on the right hand side of figure |]. They are obviously 
mapped onto each other by a rotation about the 7 -axis by 8ttGM 1 in clockwise direction. Sim- 
ilarly, r/2- is mapped onto i] 2+ by a rotation about the 7 -axis by 8nGM 2 in counter clockwise 
direction. It is shown in the appendix that the group elements representing these rotations are 

ht = e 47rGMl ^°, h 2 = e - 4nGM ^°. (2.13) 

To be precise, there are more general Lorentz rotations mapping the given edges onto each other. 
We may, for example, multiply hi from the left with a boost whose axis is parallel to the edge 
r\x-, or from the right with boost whose axis is parallel to the edge r]i + . However, there is 
yet another consistency condition. The two half spaces have to fit together as two halves of a 
single screw surface, without a kink, and this is only the case if the transition functions are pure 
rotations. 

Regarding the transition function g, which specifies how the edges A + and A_ are glued to- 
gether, there is no such additional restriction. The only consistency condition is that the Lorentz 
rotation g maps the Minkowski vector z_ onto the vector z + , hence 

z + = g x z g z =gz + g l . (2.14) 

There is a one parameter family of such group elements, and therefore there is one additional 
degree of freedom that has to be fixed. As an ansatz, we expand the matrix g in terms of the unit 
and gamma matrices, using 70, 7(0) and *y'((fr) as an orthonormal basis, 

g = ul + wjo + q~f{4>) + sj'((f)). (2.15) 



This basis is adapted to the definition Q2.1 1| ) of the vectors z±, so that the equations ( |2.14| ) 



for the coefficients u, w, q, s become as simple as possible. Additionally, we have the equation 
det(g) = 1. The resulting system of quadratic equations can be easily solved. There is a one- 
parameter family of solutions, which is given by w = and 



R 2 + {AvGS) 2 AttGQ AtcGS „ ^ 

u = -, =, q = . =, s = (2.16) 

lR 2 - (4ttGQ) 2 JR 2 - (AttGQ) 2 JR 2 - (AtcGQ) 2 



The reason for choosing the parameter Q in this particular way will become clear in a moment. 
The motivation is to have a certain symmetry between the spatial components q and s of the 
matrix g on one side, and the variables Q and S on the other side. 

Due to the square roots in the denominator, there is obviously a restricted range —R < 
A-kGQ < R for the new parameter. To avoid such a non-trivial restriction on the range of the 
phase space variables later on, it is useful to replace the variable R by a new variable X, so that 

R 2 = X 2 + (AnGQ) 2 . (2.17) 

The expression for the transition function g then simplifies slightly and becomes 

s = ffltM «*W! ( (2 , 8) 

X 
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where U > is defined so that the gr has a unit determinant, 



f/ 2 = l + (47rG) 2 ^^. (2.19) 



Here we should note the similarity to the formula fll.20p , defining the momentum vectors of the 



free particles. If the group element g is considered as the spatial momentum of the particles, the 
Q and S are generalized versions of the radial and angular momentum. 

To see what kind of group element g it is, let us rewrite it in the following alternative way. We 
replace the Cartesian coordinates Q and S by polar coordinates K and (3, so that 

AttGQ = X sinh(4vrGK) cos(/3 - 0) , AnGS = X smh(4irGK) sin(/3 - 0) . (2.20) 

This is obviously a generalization of fll.l9p , to which it reduces in the limit G — > 0. It follows 
that U = co?h(ATtGK), and the transition function becomes 

g = e^ GK ~<W = cosh(47rGir) 1 + sinh(4vrG'K) *y(j3). (2.21) 

It represents a boost with rapidity AitGK and angular direction (3. There is thus a certain relation 
between the transition functions h k and the energies M k of the particles on one side, and the 
transition function g and the spatial momentum K of the particles on the other side. Of course, 
we thereby assume that K still defines the spatial momentum, in the same way as previously for 
the free particles. 



Mass shell constraints 

To see that Mk are in fact the energies, and K is the spatial momentum of the particles, let us 
now have a closer look at the conical singularities, which arise at the corners of the half spaces 
in figure |] when they are glued together. We somehow have to ensure that the spacetime metric 
in the neighbourhood of a world line represents the gravitational field of a point particle with a 
fixed mass m^. It has to be a cone with a deficit angle of SnGrrik, when this angle is measured 
in the rest frame of the particle. In the static case, we can read off this deficit angle directly from 
the geometry of space in figure ^, or figure 0(b). In general, however, we first have to transform 
to the rest frame of the particle, as otherwise the deficit angle is proportional to the energy and 
not to the rest mass. 

An alternative way to derive the deficit angle is to consider the holonomy of the particle. 
This is, by definition, the Lorentz rotation that acts on a spacetime vector which is transported 
once around the particle, say, in clockwise direction. If the deficit angle of the particle ix k is 
SnGrrik, then the holonomy is a rotation by 8nGmk about some timelike axis. The axis is 
parallel to the world line, and thus the holonomy also defines the direction of motion of the 
particle. Massless particles can also be included. The holonomy is then a null rotation, and the 
world line is lightlike. The holonomy is in this sense a generalized, group valued momentum of 
the particle [0JJ]. 

If we think of the embedding Minkowski space in figure |as a spacetime atlas consisting 
of two charts, then a spacetime vector is simply represented by a Minkowski vector, which is 
attached to some point on one of the two half spaces. As long as we stick to this half space, the 
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parallel transport is trivial. But whenever we pass across one of the edges i] k or A from A_ to 
A + , we have to act on the vector with the appropriate transition function, hence h k or g. And 
when we pass from A + to A_, we have to act on the vector with the inverse transition function. 
There are then all together four holonomies that can be defined. We can choose the particle n k to 
be surrounded, and we can choose the half space A± in which the path begins and ends. 

Let us define u k ± E SL(2) to be the Lorentz rotation acting on a vector which is first defined 
in the half space A±, and then transported in clockwise direction around the particle n k . Sorting 
out the factor ordering and the signs, we find that 

ui + = g~ l ht, ui =h 1 g~ 1 , u 2+ = hf 1 g, u 2 =gh 2 ~ 1 . (2.22) 

For u k ± to be a rotation by 8nGm k in clockwise direction, it has to be an element of the conju- 
gacy class of e 47rGmfe " 1 ' in SL(2). This special element represents a clockwise rotation by 87iGm k 
about the 70-axis. All others are obtained by acting on this with a proper Lorentz rotation, hence 
a conjugation with some other element of the Lorentz group. The conjugacy class is defined by 
the following mass shell and positive energy condition (TT3J, [TOp 

u k = §TrK ± ) = cos(4vrGm fe ), ±Tr(u fc± 7°) > 0. (2.23) 

Since the trace of the holonomy is independent of the factor ordering, this is one equation to be 
satisfied for each particle. Actually, each particle has only one holonomy u k . The two different 
values u k+ and u k _ arise because the same physical object is represented in two coordinate 
charts. 

If the holonomy is interpreted as a generalized momentum, then the mass shell constraints 



are obviously generalizations of the free particle constraints fll.3|). There is now, however, a 
restriction on the mass parameters m k , which has no counterpart for the free particles. The 
deficit angle of a cone must be smaller then 2n, and consequently the rest mass of a particle is 
bounded from above by 1/4G = M P1 /4, where M P \ = 1/G is the Planck mass. At this upper 
bound, the cosine in ( [2.23D takes its minimum —1, and the holonomy becomes a full rotation by 
2n. For a massless particle, on the other hand, the cosine is equal to one, which means that the 
holonomy is a null rotation. 

The same restriction applies to the total energy M, which defines the total deficit angle of 
universe. We'll see later on that the energy M is in fact also bounded from below by M min = 
rrii + m 2, which is the total energy for a static state. Hence, the allowed range for the mass 
parameters is 

m k > 0, m 1 +m 2 < M P i/4. (2.24) 

To see that the free particle mass shell constraints are recovered in the limit G — > 0, let us write 
the holonomies as functions of the configuration variables. It is thereby useful to introduce the 
following notation. We define a modified set of trigonometric functions with rescaled arguments, 

fn n \ sin(27rGg) tan(27rGg) 
rag = cos(2ttG^), sn£= — ^-q — > tng = -^-q • ( 2 - 25 ) 

Some useful properties of these functions are that the relation tn g = sn g/cs g remains valid, 
and in the limit G — > we have cs g — > 1, sn g — > g, tn g —> g. We also introduce an analogous 
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set of hyperbolic functions, but this time it is useful to rescale the argument by a different factor, 



i i f a y~f \ i smh(47rG£>) tanh(47rG^) 

csn g = cosn(47rGf?J, snh g = — — , tnnp = — — . (2.26) 

47rG 4ttG 

Again, we have tnh g = snh g/csh g, and in the limit G — > we get csh g — > 1, snh £ — > g, 
tnh — > g. As an example for the application of these rescaled function, consider the definition 
( [2.20D , relating Q and S to X and (3. They simplify to 

Q = X snh X cos({3 - 0) , 5 = X snh X sin(/3 - 0) , (2.27) 

and X can be written as a function of Q, S, and X, 

snh 2 X = ^-t^. (2.28) 



According to the general rule snh g — > g, this reduces to ( |1.2ip in the limit G — > 0. Similar rules 
apply to all other formulas below. They always reduce to the free particle counterparts in the 
limit where the gravitational interaction is switched off. 

To express the holonomies and finally the mass shell constraints in terms of the energy and 
momentum variables, it is useful to make the same redefinition that we previously also made for 
the free particles. We replace the mass parameters m k by a total mass /i and a relative mass u, 
assuming without loss of generality that m 2 > mi, 

/i = m 2 + mi, v — m 2 — m ii < v < \l < Mpi/4. (2.29) 

Like for the free particles, the special cases are v = 0, where both particles have the same mass, 
and v — n, where at least one of the particles is massless. The only new feature is the upper 
bound for fi, which goes to infinity in the limit G — > 0. Similarly, we also replace the energy 
variables M k by the total energy M and a relative energy V, so that 

M = M 2 + M U V = M 2 -M U \V\ < M < M P1 /4. (2.30) 

The upper bound on M follows from the fact that the total deficit angle of the spinning cone 
must be smaller than 2n. The restriction on V is the positive energy condition, which implies 
that M k > for both particles. 

Using all this, we can finally express the transition functions ( [2.1 3| ) and Q2.21| ) in terms of the 
new energy and momentum variables. What we get is 

h x = cs(M - V) 1 + 2ttG sn(M - V) 7o , 
h 2 = cs(M + V) 1 - 2nG sn(M + V) 7o , 

g = csh XI + 4ttG snh K~/([3). (2.31) 

It is then not difficult to evaluate the traces of the holonomies, and to express the mass shell 
constraints in terms of the same variables. They are given by 

mi = cshXcs(M - V) = cs(fj, - u), u 2 = cshXcs(M + V) = cs(> + v). (2.32) 
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This is the same as (2.36) in [{131], if K is interpreted as the common spatial momentum of both 
particles in the centre of mass frame, and (M ± V)/2 are the energies of the particles, which are 
in general different. Finally, we have to rescale the actual mass shell constraints by an appropriate 
power of G, in order to get the correct limit G — » 0, 

_ cshifcs(M - V) - cs(/i - v) _ cshif cs(M + V) - cs(/x + i/) 

Expanding this up to the second order in G, we recover the free particle constraints fll.lip in the 
limit G — > 0. And the same applies to the positive energy condition ( |2.30D , where the upper 
bound on M goes to infinity. 

We can then proceed in the same way as before. We use one of the mass shell constraints to 
eliminate the relative energy V. There is then a single mass shell constraint left, which is going to 
be the generator of the time evolution. And we also get an even dimensional phase space which 
is spanned by six independent variables. The first step is again to define the linear combinations 

, csh K sn M sn V — sn /i sn v 
V = C 2 -Ci = , 

„, csh K cs M cs V — cs u cs v _ , 

£=C * + C ' = 2(2,0)' (234) 

Here we used some trigonometric identities to simplify the results. In the limit G — > 0, these 



constraints are equal to (1.12). However, it is now a little bit more complicated to decouple the 



variables K and V. We have to take another linear combination, namely 

V _V + (2vrG) 2 tnM tnV £' _ tn fi tnu - tnMtnV 
cs n cs v 2 

Note that the functions 1/cs/i and 1/cs v, and similarly tn M and tn 1/ are well defined within 
the range Q2.29D and ( |2.30| ). The first poles of these functions are just outside the allowed range 



of the arguments. Now, the positive energy condition implies tnM > 0, and the function tn is 
invertible within the allowed range of V . We can solve the equation V = for V, setting 



/tn a tn zA 

V = atn " aj . (2.36) 



V tnM 



The function atn is the inverse of tn, 



atn q = BXct ^ 2 ^ G ^> _^ -Mpi/4 < atn g < M Pi /4, (2.37) 
2nG 

and in the limit G — » we have atn g ^ g. Using this and the fact that tn and atn are both 
monotonically increasing, the positive energy condition becomes a non-trivial condition to be 
imposed on M, namely 

tn u tn v / , \ 

| tn V| = < tnM ^> atn (ytn/itn u j < M < M P1 /4. (2.38) 

This is obviously a generalized version on fll.151 ), with all mass parameters and energy variables 
replaced by their rescaled tangent. This is in fact a rather general rule, which tells us how to 
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obtain the various structures of the Kepler system from those of the free particle system. Consider 
for example the mass shell constraint that remains after eliminating V from ( |2.34D . To derive this, 
we have to take yet another linear combination, namely 

csh K cs M cs V + cs u cs v . csh K sn M sn V + sn a sn v . 

£ - 2^m £+ 2*? M v - <239) 



Again, the coefficients are well defined, because the range ( |2.30D of M is exactly the interval 



where both sn M and cs M are different from zero. Evaluating this linear combination, V drops 
out, and after some trigonometric simplifications we get 

£ = snh 2 K - A 2 F(M) « 0. (2.40) 

Here, < A < 1 is a constant, which only depends on the mass parameters fj, and u, and F is the 
same function that also appears in ( ]1.14D , but once again with the masses /i, u, and the energy M 
are replaced by tn /i, tn v, and tn M, respectively, 

, , (tn 2 M - tn 2 a) (tn 2 M - tn 2 u) 
A = cs/zcs I /, F(M) = { - 1. (2.41) 

Roughly speaking, the rule is to replace all mass and energy quantities g by tn g, and to replace 
the momentum K by snh K. Then we obtain the mass shell constraints of the Kepler system from 
those of the free particle system, and also, for example, the relation ( [2.28D between the radial 



momentum Q, the angular momentum S, and the total spatial momentum K of the particles. As 
all these modifications become identities in the limit G — > 0, we can say that the Kepler system 
is a deformation of the free particle system, and Newton's constant is the deformation parameter. 
We'll now see that this also applies to the various phase space structures, such as the symplectic 
potential, the Poisson brackets, and the Hamiltonian. 

The phase space 

So far, we have seen that all physically relevant quantities, including the geometry of space, the 
positions of the particles with respect to the centre of mass frame, as well as the holonomies and 
the mass shell constraints, can be expressed as a function of six independent variables, the energy 
and momentum variables M, Q, S, and the time and position variables T, X, <p. In analogy to 
the free particles, we define the extended phase space 

V = { (M,Q,S;T,X,<f>) | X > 0, = + 2tt }. (2.42) 

Note that we do not impose any restriction on the energy M at this point. The positive energy 
condition ( |2.38[ ) is imposed together with the constraint ( |2.40| ), and this defines the physical 



subspace. The Hamiltonian on V is again proportional to the mass shell constraint. Thus if we 
introduce a multiplier £, we have in analogy to ( |1.24| ) 

H = (£. (2.43) 

To derive the equations of motion, we also have to know the symplectic structure on V. But so 
far, there is no natural way to define it, and there is also no reason why it should be the same 
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as that of the free particle system. Since we do not want to make any additional assumptions, 
we can only derive the symplectic structure from the underlying field theory of Einstein gravity. 
Thus, we have to apply a straightforward phase space reduction to the Einstein Hilbert action. 

The actual derivation is rather involved and technical, but on the other hand it can be carried 
out without further complications for a general multi particle model. It is therefore given in a 
separate article [JTQt] . We are here not going to say anything in detail about this derivation. We 
just take the general result, and adapt it to the special case of a two particle system. The general 
expression (3.9) in JlOt], adapted to the special situation given in figure |], is 

6 = - £ M k (dT k + AGS d<p k ) - -L Tr^-Mf? z+). (2.44) 

k o7TCt 

Here, M k are the energies of the particles ix k , T k and <p k are their absolute positions with respect 
to the conical reference frame, z + is the relative position vector of the particles in the half space 
A + , and g is the transition function mapping the edge A_ onto A+. 

This expression does not look symmetric, but in fact it is, which can be seen by replacing z + 
and g with z_ and g 1 . To express the symplectic potential in terms of our six independent phase 
space variables, we have to insert the expressions ( [2.10D for T k and 4> k , ( |2.1 1| ) for z±, and ( |2.18| ) 
for g. The result is very simple and reads, up to a total derivative that can be neglected, 

= X~ l Q dX + (1 — AGM) Sd(f) — M dT. (2.45) 

This is almost the same as the free particle expression fll.22p . There is only one modification. 
The angular momentum S, which is conjugate to the orientation of the particles, is rescaled 
by a factor 1 — AGM. This factor takes values between zero and one, and obviously it has to 
do with the conical geometry of the spacetime at spatial infinity. As a consequence, we have the 
following almost canonical Poisson brackets 

{M,T} = -1, {Q,X} = X, (2.46) 

but the following non-canonical brackets involving the angular momentum S, 

1 AG 9 

^ = T^GM' {S - T} = -T^GM- (2 ' 47) 

The last one results from an off-diagonal term in the symplectic two-form Q = d0, involving a 
product of dM and dep. The brackets can most easily be derived if we notice that the canonically 
conjugate angular momentum is actually 

J = (1 - AGM) S {J, <f)} = 1. (2.48) 



All other brackets with J are zero. The brackets ( [2.47D are then easily found by expressing S 



as a function of J and M. But let us nevertheless stick to S as one of the basic phase space 
variables, because it has an immediate geometric interpretation, defining the time offset of the 
spinning cone. It also shows up in the transition function g defined in ( |2.18| ), and in the mass 
shell constraint ( [2.40D , implicitly through the definition ( [2.28D of K. 



Now we have all the structures at hand which we need to describe the kinematical and dy- 
namical features of the Kepler system. In particular, we can now derive the equations of motion, 
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study the classical trajectories, and finally we can quantize the Kepler system. At each step, the 
result will be a deformed version of the corresponding free particle result. The most interesting 
feature is thereby the deformation of the symplectic structure by the factor in front of the angular 
momentum. Apparently, this is only a marginal modification. However, it turns out that it is 
responsible for some unexpected effects at the quantum level. For example, the particles can no 
longer be localized in space, and they cannot come closer to each other than a specific distance. 

But for the moment we shall stick to the classical phase space. Some basic features of the 
Kepler system can be inferred immediately from the given phase space structures. For example, 
we have a two dimensional rigid symmetry group of time translations T i— > T — AT and spatial 
rotations <p |— ► + A0. The Hamiltonian and the symplectic potential are both invariant. At the 
spacetime level, these are the Killing symmetries of the spinning cone. The rigid symmetries are 
the possible translations and rotations of the universe with respect to the reference frame [[14]]. 
The conserved charges associated with these symmetries are obviously M and J. 

Actually, we expected S to be the angular momentum. It is a function of M and J, thus also a 
conserved charge, but the associated symmetry is a combination of a time translation and a spatial 
rotation. It is a screw rotation of the spinning cone. The Killing vector of this symmetry is a linear 
combination of the Killing vectors of time translations and spatial rotations. It is orthogonal to 
the Killing vector of time translations, whereas the rotational Killing vector associated with J 
has closed orbits of affine length 2n. Unless the spinning cone is static, this is not the same. The 
definition of an angular momentum is therefore somewhat ambiguous. 



Complete reduction 

We can also define a reduced phase space V, in the very same way as before. Instead of a mass 
shell constraint, the time evolution is then provided by a Hamiltonian that represents the physical 
energy. To go over from the constrained to the unconstrained formulation, we have to solve 
the mass shell constraint £ = 0, and additionally we have to impose a suitable gauge condition. 
As all this is completely analogous to the free particles. We can choose the same natural gauge 
condition, requiring the ADM time t to be equal to the absolute time T in the centre of mass 
frame. What we get is a pair of second class constraints, 

snh 2 K = A 2 F(M), T — t. (2.49) 



Since F is more or less the same function as before, just with some modified mass parameters, it 



can still be inverted, defining M as a function of K. The explicit solution is given in ( [2.53D below. 
It is thus possible to remove the canonical pair (M, T), and to go over to a four dimensional, 
completely reduced phase space 



V = { (Q,S;X,, 



X > 0, 



2tt }. 



(2.50) 



To derive the symplectic potential and the Hamiltonian H on V, we have to start from the 
extended symplectic potential on V, again because the gauge condition above is explicitly time 
dependent. According to ( |2.43[ ) and ( |2.45| ), it is given by 



G - Hdt = X^QdX + (1 - AGM) S d(f> - MdT - (£ dt. 



(2.51) 
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On the reduced phase space, this becomes 



6-Hdt = X~ l Q dX + (1 - AGH) S d0 - H dt. (2.52) 



To derive the reduced Hamiltonian H, we have to solve the equation Q2.49D for M, so that the 
positive energy condition (|2.38|) is satisfied. The unique solution is 



H = F~ l [A- 2 snh 2 K) = atn ( J A~' 2 snh 2 K + m 2 + J A~ 2 snli 2 K + m 2 2 ) , (2.53) 



where K is given by ( ]2.28| ) as a function of the Q, S, and X. The deformed mass parameters rrik 



are just some useful abbreviations. They are given by 

tnu — tnz/ „ tnu + tnz/ 
mi = — ^ , m 2 = — ^ . (2.54) 

So, we find a somewhat deformed Hamiltonian, as compared to the free particle definition ( |1.31| ). 
It still depends only on the momentum K, and not, for example, on the relative position X. This 
we might have expected as a kind of gravitational potential. However, in three dimensional Ein- 
stein gravity, there are no local gravitational forces, and therefore we do not have a gravitational 
potential. But nevertheless, there is something like a gravitational binding energy. To see this, 
let us briefly discuss the range of H, and the way it depends on the momentum K. First of all, 
we easily see that H is minimal for K = 0, where it takes the value 

M min = mi + m 2 - (2.55) 

As for the free particles, the static states with K = are excluded if at least one of the particles 
is massless. In this case, the positive energy condition 02.38D requires M > M min . Moreover, 



for small momenta K, the energy also increases linearly if a massless particle is present, and 
quadratically if both particles are massive. For massive particles, we can also consider a non- 
relativistic limit, where K is small and the particles are moving slowly compared to the speed of 
light. Expanding ( |2.53j ) up to the second order in K gives 



1 ., ~ cs 2 v tn 2 u — sn 2 v 

M w M min + — K 2 , with m = — . (2.56) 

2m 4tn/i 

Apart from the deformed expression for the reduced mass m, we have the usual non-relativistic 
behaviour of the kinetic energy. It can be shown that m > m, thus the reduced mass of the 



coupled system is always larger than the reduced mass of the free particles in ( ]1 .33| ). This is the 
only effect of the gravitational coupling for low momenta. For a given value of K, the energy 
of the coupled system is slightly smaller than the free energy. There is still no gravitational 
potential. In particular, we do not recover Newtonian gravity. In three spacetime dimensions, 
this is not the non-relativistic limit of Einstein gravity. The relation between the two theories 
only exists in higher dimensions. 

The behaviour of H for large momenta K is changed more drastically, as compared to the 
free particle system. Asymptotically, the energy is no longer linear in K. The right hand side 



in fl2.53| ) is bounded from above by Mpi/4. The energy asymptotically approaches the upper 



bound for large momenta K. This is also a non-trivial result, because we might have expected 
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Figure 5: The energy momentum relation for the Kepler system (solid lines), and the free 
particles (broken lines). For low momenta, the behaviour is almost the same. For large mo- 
menta, the energy of the free particle system increases linearly with the momentum, whereas 
that of the coupled system approaches the Planck energy from below. All quantities are given 
in Planck units Mpi = l/G, and the upper bound for the total energy is M max = Mpi/4. 



that the upper bound Mpi/4 for the total energy M imposes some upper bound on the spatial 
momentum K as well. But it seems that the gravitational interaction knows about this upper 
bound, and takes it into account in the relation between momentum and energy. To compare the 
energy momentum relation for the free and the coupled particles, we have plotted them for some 
typical mass parameters in figure |5|. 

For low momenta, the behaviour is very similar, we just have a slightly modified reduced mass 
m > m for massive particles. This implies that the energy of the coupled system is a little bit 
smaller than the free particle energy. For large momenta, however, we see that the energy of the 
coupled system approaches M max = Mpi/4, whereas the energy of the free particles increases 
unboundedly. We can think of the difference between the free particle energy and the total 
energy of the coupled system as the negative energy of the gravitational field. This obviously 
compensates for the increasing kinetic energy of the particles, so that the total energy remains 
below the upper bound, which a quarter of the Planck energy. 

If the momentum K is of the order of one tenth of the Planck scale, then the negative binding 
energy is already of the same order of magnitude as the free particle kinetic energy. For momenta 
which are larger than the Planck scale, the energy has more or less saturated at the upper bound. 
The same behaviour of the energy can be found for a single particle system [Op, where the 
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upper bound is Mpi/8, thus half of the maximal energy here. This is a generic feature of three 
dimensional general relativity. It has to do with its topological nature. The gravitational field 
does not fall off at spatial infinity, which means that the spacetime is not asymptotically flat. 
Instead, it is asymptotically conical, and therefore it is not possible for the universe to contain 
more energy than M max = Mpi/4. 

But now, let us turn to the more interesting object on the reduced phase space V, the symplectic 



structure. It follows from ( [2.52D that the reduced symplectic potential is given by 

6 = X~ 1 QdX + (1 -AGH)Sd(j). (2.57) 

The somewhat curious feature of this expression is the Hamiltonian H, which appears at the place 
where originally the energy M has shown up. Since H depends on Q, S and X, the resulting 
Poisson brackets become very complicated, because the two form i? = d<9 is far from being 
diagonal in the given variables. 

It is therefore more appropriate first to look for a set of canonically conjugate variables and 
then to define the Poisson brackets. Clearly, what we have to do is to replace the variable S by 
the canonical angular momentum J = (1 — AGH)S. Doing so, the symplectic potential and the 
resulting Poisson brackets become 

6 = X^QdX + Jd<f) {Q,X} = X, {J,<f)} = l. (2.58) 

All other brackets of Q, J, X, and 4> are zero. We can further replace X and Q by the geodesic 
distance R of the particles, given by ( |2.17| ), and its canonically conjugate momentum P. The 
transformation reads 

R 2 = X 2 + (AnGQ) 2 , snhP = %, (2.59) 



X 



and the inverse transformation is 



X = ——, Q = Rtnh P. (2.60) 
csii P 

Inserting this into the symplectic potential gives 

6 = PdR + Jd(f) {P,R} = 1, {J,(f)} = l. (2.61) 

The chart (P, J; R, <fi) is thus canonical on V, and the coordinates have a straightforward physical 
interpretation. This is also the canonical chart introduced in pOQ. 

In the chart (Q, J; X, <p) the Hamiltonian H is determined implicitly as the solution to 

A 2 F(H) = ^ K = 1^1 = gl+O^OT , (Z62) 



and the corresponding equation in the chart (R, 0, P, J) follows by the substitution ( |2.6(jp . As 
F(H) involves the rescaled tangent of H, the equation is transcendental in either chart. The price 
for simplifying the symplectic potential fl2.57p into ( |2.58| ) or into the fully canonical form ( |2.61| ) 



therefore is that H cannot be expressed as an elementary function of the variables. 
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Trajectories 



To find the classical trajectories, we return to the constrained formulation based on the six di- 
mensional phase space V, with the additional independent variables M and T, the symplectic 
structure ( |2.45| ), and the Hamiltonian constraint ( [2.43[ ). As V is the gravitating version of the 
special-relativistic phase space of section [j], we can follow the analysis of section [j] with the 
appropriate deformations. 

The starting point is the Hamiltonian, which generates the time evolution with respect to the 
unphysical ADM time t. It is given by the mass shell constraint ( [2.40D , and a multiplier £, which 
is some arbitrary function of t, 

H = (£, £= Q2 ^f -A 2 F(M). (2.63) 

This is almost identical to ( |1.14{ ), and we can now use the brackets Q2.46| ) to derive the time 
evolution equations. The total energy M and the angular momentum S are of course conserved 
charges, and the same holds for the total spatial momentum K, 

M = {H, M} = 0, S={H 1 S} = K = {H, K} = 0. (2.64) 

The time evolution equations for the radial coordinates are almost the same as ( |1.35[ ), 

X = {H,X} = 2(%, Q = {H,Q} = 2( Q2 + S2 =2Csnh 2 K. (2.65) 



The only essential modifications arise when we derive the evolution equations for <fi and T. They 
differ from fll.36| ), because of the non-canonical brackets ( |2.47| ) of S with <p and T. For the 
angular orientation, we find 



^-^ (i-Li)x^ (2 - 66) 



and for the clock T we get 



± = {H,T) = C,(^nM)- T ^M^. (2.67) 

At first sight, these equations of motion are somewhat complicated, but it is still possible to solve 
them explicitly, in the same way as in section [T|. First, we observe that M, S and K are constants 
of motion. Then, we introduce a function e(t) so that e(t) = 2((t). Using the same arguments 
that led to Q1.37[ ), we conclude that 

M{t) = M , S(t) = So, K(t) = K Q , Q{t) = e(t) snh 2 K , (2.68) 

where M , S and K are integration constants. It is then straightforward to solve the equation 
of motion for X. The solution is 



X(t) = JR 2 + t\t) snh 2 K , (2.69) 
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where R is some integration constant. And yet another constant <fi arises when we then solve 
the equation of motion for </>, 

. . 1 ( / \ snh 2 K \ 

^ ) = ^ + r34G^ arctan ( E(t) ^j- (2J0) 

And finally, we can also solve the equation of motion for T. The solution is 

T(t) = T + e(t) T -^ Wo arctan [e(t) j , (2.71) 

where T is yet another integration constant. All together, the resulting trajectory is parameter- 
ized by the gauge function e(t), and six integration constants M , K , S , T , R , (fi . Again, 
only four of them are independent. It follows from the constraint S = 0, and the definition of K 
as a function of Q, S, and X that 

snh 2 K = A 2 F(M ), R snhK = \S \. (2.72) 

These relations and all the formulas above reduce to their free particle counterparts in the limit 
G — > 0. What applies to the phase space structure, also applies to the classical trajectories of the 
Kepler system. They reduce to those of the free particles when the gravitational interaction is 
switched off. This is also a non-trivial statement. It implies, for example, that there are no bound 
states, again as a consequence of the absence of local gravitational forces. 

As seen in section [I], the trajectories on the reduced phase space V are obtained from those 
on V by imposing the gauge condition T{t) = t. We shall now verify that this gauge is always 
accessible and unique. By construction, the trajectories on V then solve Hamilton's equations 
with the Hamiltonian defined implictly by ( |2.62[ ). We have to show is that the right hand side of 



Q2.71D is a monotonically increasing function of e(t), and that its range is the whole real line. If 
this is the case, then the equation T(t) = t can always be solved for e(t). 

That the range is the whole real line can be seen quite easily. In the range of M which is 
allowed by the positive energy condition, F'(Mq) is positive. Hence, if e(i) is large, then the 
arctan term can be neglected as compared to the linear term, which is unbounded. So, what 



remains to be shown is that the derivative of the right hand side of ( |2.71| ) with respect to e{t) is 



positive. We know already that, up to a factor of 2£ = e, this derivative is equal to the right hand 



side of ( |2.67| ). This is how it was derived. Hence, we have to show that for all physical states 

8G S 2 
1 - AGM ~X r 



or* c2 

^ 2 F'(M)>-^-" (2.73) 



The constraint ( [2.63D implies that 



X 2 ~ X 
It is therefore sufficient to show that 



< g Z_g = A 2 F(M). (2.74) 



F ' (M)> r^M F(M) - (2 - 75) 
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This can be simplified further, using the following estimation. For < a < tt we have tt — a > 
sin a. With a = AttGM, this implies 

l-4 G M> Sin < faGM >=4 G snMcsM - g ° < » -. (2.76) 

tt 1-4GM snMcsM 



Hence, it is also sufficient to show 



F'(M) > J^L ,2.77) 
sn M cs M 



Now, all we have to do is to insert the definition <J2.41|) for F and compute its derivative. It is 



then very easy to verify that this inequality is indeed satisfied, provided that the positive energy 
condition holds, that is tn M > y/tnfitn 7J. Thus, we conclude that a gauge function e(i) exists, 
so that T(t) = t, although we do not know it explicitly. 

Scattering 

Let us now describe the physical properties of the trajectories, in particular the way they deviate 
from the free particle trajectories. Since the relative position X is still an auxiliary variable, we 
should first replace it by the actual geodesic distance R between the particles. And it is also 



useful to consider the previously defined canonically conjugate momentum P. Using Q2.59D , we 
find that on the given trajectories we have 

e(t) snh 2 K 



R(t) = ^R 2 + le 2 (t)snh\2K ), snhP(t) = w u (2.78) 

'R 2 + e 2 (t)snh 2 K 



Here, we can first of all see that the assumption was correct that K represents the ingoing and 
outgoing momentum of the particles. If R(t) is large, then the canonically conjugate momentum 
P(t) approaches ±K . As for the free particles, we expect the eigenvalues of K to specify the 
inverse radial wavelength at infinity, when the system is quantized. 

But now, let us consider the actual trajectories. First we assume that So ^ 0. According to 
( |2.72| ), this implies that R > and K > 0. In this case, the particles approach each other 
from infinity from some direction in , they reach a minimal distance R at the absolute time 
T and with angular orientation O , and afterwards the separate again, moving to infinity into a 
direction out . The interpretation of the integration constants M and So is also obvious. These 
are the total energy and the angular momentum, which are constant along the trajectory. So far, 
everything is the same as before. 

However, there is now a crucial difference, because now we have a real scattering. The angular 
directions in and out are given by 

It is no longer so that out = 0; n ± n. Instead, the scattering angle is 

, , . sgn So 
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The scattering angle depends on the deficit angle of the spinning cone, thus on the energy of the 
particles, and the direction of the scattering depends on the sign of the angular momentum, thus 
on which side the two particles pass each other. This is not a surprise, because it is the typical 
behaviour of a pair geodesies on a cone, and of course also on a spinning cone. 

In fact, the world lines of the particles can effectively be considered as a pair of geodesies on 
antipodal sides of the spinning cone [EO]. However, the crucial point is that the geometry of this 
spinning cone depends on the energy and the angular momentum of the particles. The particles 
are not just moving on a fixed background cone. It is also interesting to note that for large 
momenta K , thus for M ~ M max , the factor in the denominator in ( |2.80[ ) becomes arbitrarily 
small. In this case, the particles wind around each other many times before they separate again. 
We also see that the scattering angle is always larger than n. This means, in a sense, that gravity 
is attractive, even though there are no local forces. 

The case S = is again somewhat special. According to ( |2.72[ ), we either have K = 
or R = 0. For K = 0, we have a static state. The particles are at rest, with a fixed distance 
Ro > 0, and a fixed angular orientation <fr . The geometry of space is that of figure [|, and it is 
constant in time. The total energy is M = M min , and the integration constant T is redundant. 
The more interesting case is So = and R = 0. Let us once again consider the limit So — > 
and Ro — * 0, with K > and all other integration constant fixed. What we find is similar to 

Cgg ), 

m = m\ s -^, m = 00 + f sg ne( t). (2.8D 

Again, the particles approach each other on a radial line with direction in , then they hit in a head 
on collision at the absolute time T , and after that they continue on a radial line with direction 
4> out . However, there is now a crucial difference to the free particles. The relation between </> in 
and (fi out is no longer independent of the sign of S , and thus the limit is actually not well defined. 
We have 

, M J 0in for T < T , vr 

0(t) = U out forT>T , ^^ITlGMo' (2 " 82) 

which depends on the sign, unless 1/(1 — 4GM ) £ Z. These are very special cases where the 
deficit angle of the universe takes a value in a certain discrete set. 

Clearly, this is the typical behaviour of a geodesic on a cone. If the geodesic hits the tip, then it 
does not have a unique continuation. If it is considered as a limit of a family of geodesies missing 
the tip, then its continuation depends on the way we take the limit. To get a unique time evolution, 
we have to say explicitly what the particles have to do when they hit each other. Quite remarkably, 
this ambiguity disappears at the quantum level. The following semi-classical argument explains 
why. The problem only arises when the angular momentum S vanishes. But then, the angular 
orientation <fi is smeared out uniformly over all directions. Hence, for 5 = there is anyway no 
well defined scattering angle, because both in and 4> out are maximally uncertain. 

Apart from this ambiguity, the classical dynamics of the Kepler system is well defined on the 
extended phase space V, and the trajectories are easy to derive. Implicitly, this also defines the 
trajectories on the reduced phase space V. We know that they exist, and of course they have the 
same physical properties, but we cannot write them down explicitly. To go back to the spacetime 
picture, we have to perform the construction in figure |] at each moment of ADM time t. What 
we get is a space that evolves in time, and this provides a foliation of the spacetime manifold 
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by slices of constant absolute time in the centre of mass frame. We can either take this general 
relativistic point of view, but we may also stick to a simplified point of view where we just think 
of two particles moving in a two dimensional space. 



3 The quantum Kepler system 

With all the preparations made in the previous sections, we can now, finally, quantize the Kepler 
system. As an ansatz, we shall first try the Schrodinger method, based on the reduced phase 
space V. This turns out to fail, due to a somewhat peculiar feature of the classical phase space. 
Its global structure is not as simple as it appears, which implies that a straightforward choice of a 
set of commuting configuration variables does not provide a proper basis of the quantum Hilbert 
space. We already noticed that the definition of the Poisson brackets and the Hamiltonian was 
technically somewhat involved, and we were not able, for example, to define the Hamiltonian 
explicitly as a functions of a set of canonically conjugate variables. 

This technical problem turns into a principle one at the quantum level. The way out is to 
use the Dirac quantization method, based on the extended phase space V. On this phase space, 
the mass shell constraint can easily be expressed as a function of a set of canonically conjugate 
variables. Instead of solving the Schrodinger equation, we have to impose a constraint. It be- 
comes a generalized Klein Gordon equation and defines the physical Hilbert space. Just like the 
classical phase space, the quantum Hilbert space turns out to be a deformed version of the free 
particle Hilbert space. Finally, we are going to write down the energy eigenstates explicitly as 
wave functions in space and time, and from them we read off some interesting effects of quantum 
gravity. 



Schrodinger approach 

The most straightforward way to quantize the Kepler system seems to be a Schrodinger quan- 
tization, based on the completely reduced phase space V defined in ( |2.50j ), and spanned by the 



variables Q, S, X, and <f). Apparently, the only technical problem is to define the Hamiltonian 
H explicitly as a quantum operator, which is given by ( |2.53[ ). However, a first obstacle already 
arises when we look at the symplectic potential fl2.57| ), which explicitly depends on this function, 

= X~ 1 QdX + (1 -AGH)Sd(j). (3.1) 

This we have to take into account when we set up the operator representation. We saw, however, 
that we may equally well replace the variable S by 

J = (1 - AGH) S 6 = X~ 1 QdX + Jd0, (3.2) 

and consequently we have the canonical brackets 

{Q,X}=X, {J,0} = 1. (3.3) 

Now the quantization is straightforward. We choose a basis where X and <p are diagonal. The 
wave functions is ip(x, ip), with x > 0, and ip has a period of 2n if the particles are distinguish- 
able, or the period is n if the particles are identical. We introduce a statistics parameter A, in the 
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same way as previously for the free particles. We may also define the same scalar product, 



(^1^2)= / xdxdip^x, if) ip 2 (x, p), (3.4) 



which implies that ip(x, cp) is the usual probability amplitude in polar coordinates. And it is also 
straightforward to define the operator representation, which is given by 

Xip(x,ip) = xip(x,(p), <j>i()(x,<p) = (pil)(x,(p), (3.5) 

for the position operators, and the momentum operators are 

d - d 

Q ip(x, <f) = -ih — xip(x, <p), J ip(x, <p) = -ih — if>(x, <p). (3.6) 

So far, everything is the same as before. It seems that the only technical problem is to express 
the Hamiltonian H in terms of X, Q, and J. It is now implicitly defined by the equation Q2.62| ), 
thus 

Q 2 + (1 - AGH)- 2 J 2 



X 2 



A F(H). (3.7) 



Since we could not even solve this equation explicitly at the classical level, it is of course also 
impossible to derive an explicit operator representation for H at the quantum level. Nevertheless, 
since we know that the Hamiltonian is well defined on the classical phase space, this only seems 
to be a technical problem. But if we look at all this more closely, we find that the problem is 
actually not a technical, but a principle one. In fact, we made a wrong assumption in the above 
derivation. 

To set up a quantum representation like the one given above, it is not sufficient to pick out a 
complete set of commuting phase space variables like X and <p, and to define the wave function 
ip(x, p) to be a function on this configuration space. There is a second consistency condition, 
which is sometimes overlooked because it is usually immediately obvious. The classical phase 
space has to be the cotangent bundle of the configuration space, equipped with its canonical 
symplectic structure. At each point (X, <p) of the configuration space, the canonically conjugate 
momenta, which are (Q/X, J) in our case, have to span a vector space. Splitting the phase space 
into a configuration and a momentum space like this is called a polarization [[Z3[. 

Without going into any details, let us roughly show what goes wrong when the second criterion 
is not satisfied. Given a Hilbert space such that to each point x in the configuration space there 
corresponds a basis state ijj x , then for every bounded function F on the configuration space there 
exists a bounded self adjoint operator, whose eigenstates are ip x with eigenvalues F(x). This is 
just usual quantum mechanics, for example in the position representation. Consequently, there is 
also a family of unitary operators e lrF ^ h , representing the flow of F. On the classical phase space, 
this flow is a translation in the conjugate momentum space. Therefore, the whole construction is 
consistent only if the momentum space at each point x of the configuration space is at least an 
affine space. It can be made a vector space by choosing some origin. 

Now, coming back to the Kepler system, let us check whether the second criterion is satisfied 
or not. We have to find out whether for every fixed point in the configuration space (X, 0), the 
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canonically conjugate momenta (Q/X, J) span a vector space. Unfortunately, this is not the 
case. To see this, consider the following estimation, 



^ = (1 - AGE) 2 ^ < (1 - AGH) 2 S2 + 2 Q2 = (1 - 4G7J) 2 A 2 F(H). (3.8) 
The first factor can be estimated further by 

1-AGH< 1 (3.9) 
7r 2 G tn/i 

which follows from cot a > (7r/2 — a) for < a < 7r/2, and with a = 2nGH. Using this, we 
find that 

J 2 A 2 F(H) A 2 

— < — — < (3 10) 

X 2 vr 4 G 2 tn 2 # 4vr 4 G 2 ' ^ ; 



The last inequality follows from the following property of the function F defined in Q2.41| ). For 
tn H > tn /i > tn u, which follows from the positive energy condition, we have 

F(H) _ (tn 2 H - tn 2 /i) (tn 2 H - tn 2 v) 1 

We~ ax^Fh < 4' ( } 

Hence, all together it follows that 

I-/! < Jma*, where J max = 7 ^ 7 ;X. (3.12) 

The canonical angular momentum J is restricted to some interval around zero, whose size is 
determined by the configuration variable X. Thus, at a fixed point (X, 0) of the configuration 
space, the canonical momenta (Q/X, J) do not span a vector space. We have to conclude that the 
above definition of a position representation is not the correct quantization of the given classical 
phase space. 

Of course, we could now argue that we simply picked out the wrong configuration space. We 
just have to find an appropriate configuration space to perform the correct quantization. But then 
we would lose the simple physical interpretation of the configuration variables, and consequently 
the direct interpretation of the wave function as a probability amplitude in the usual sense. There- 
fore, to avoid this unnecessary complication, we shall not try to quantize the Kepler system bases 
on the reduced phase space V 

Dirac approach 

The whole problem can be circumvented when we use the Dirac method instead. The extended 



phase space V defined in (2.42) has the additional independent variables M and T. It is a proper 



cotangent bundle with its canonical symplectic potential ( [2.45D , 

= X~ X Q dX + (1 - AGM) Sd(p - MdT. (3.13) 

The configuration space is spanned by (T, X, <p), and the canonically conjugate momenta 
(— M, Q/X, (1 — AGM)S) have an unrestricted range for fixed T, X, and 0. The quantization is 
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straightforward. We have a wave function tp(r, x, cp), with x > 0, and regarding the periodicity 
in cp we have the same relations as before. If the particles are distinguishable, we have 

iP(r,x,<p + 2ii) =e 27TiX 4>(r,x,<p), (3.14) 

and for identical particles the stronger relation is 

^(r.x^ + n) =e 7TiX i/j(r,x,i P ), (3.15) 

where the statistics parameter A is a fixed real number. The position operators are given by 

Xip(r,x,(p) = xip(r,x,(p), $ip(T,x,<p) = cpip(r,x,(p), (3.16) 

where the last equation is again to be understood with </> replaced by a periodic function thereof. 
For the clock and the canonically conjugate energy we have 

^ d 
Ti/i(T,x } (p) = Tifj{T ) x,(p), Mi/;(t,x,<p) = ih — ij(r,x,ip). (3.17) 

And finally, we need the momentum operators Q and S. As S can be expressed as a simple 
function of M and J, let us first define 

d ~ d 

Qil>(T,x,(p) = -ih — xi/;(T,x,<p), J4>(T,x,ip) = -ih — 4>(T,x,p). (3.18) 

The operator for S can then be expressed as 

Si(>(T,x,tp) = -m (l-Ai£^j ^-1>(r,x,<p). (3.19) 

The constant £ = GU is the Planck length, which is here actually interpreted as a Planck time. 
So, we have a somewhat unusual operator representation for the angular momentum S, which is 
non-local in r. But this is not a serious problem. Finally, we have the following scalar product, 
with respect to which all these operators are self adjoint, and which is formally the same as 

( TO 

(■01 1^2) = J dr xdxdipi>i(T, x, ip) 2 (t, x, <p). (3.20) 

This completes the definition of the extended Hilbert space. Now we have to impose the mass 
shell constraint to define the physical Hilbert space. The constraint is given by ( |2.63[ ), so we 
should first try to diagonalize the operator 

snh 2 K = X- 1 (Q 2 + S 2 ) X' 1 . (3.21) 

It acts on a wave function as 

snh 2 K ip(r, x, if) = —h 2 A 0(r, x, if), (3.22) 
where A is a deformed Laplacian, 

A = x -*± x ± + ^(l-m£fJL±. (3.23) 

OX OX \ OT J Oif Oif 
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Now, suppose that the operator 1 — 4GM in parenthesis is also diagonal. Then A is the Laplacian 
on a cone with a deficit angle of 8tcGM. This is not surprising. We found that the classical parti- 
cles are effectively moving on such a cone. The eigenfunctions are still Bessel function, however 
with slightly modified indices and arguments. The normalized eigenstates are parameterized by 
three quantum numbers co, k, and s, and they are explicitly given by 



*<*,M;r,,,,) = ^e-e^ . (3,4) 



This reduces to ( |1.67D in the limit G — > 0, which implies £ — * and snh g — > g. We are here 



imposing the same regularity condition at x = that we also imposed in section [I]. Following the 
arguments in [|22|], we require the wave function to stay finite. This makes ( |3.22[ ) a well defined 



self adjoint operator. The eigenstates are normalized in the same way as ( J1.68D , thus 

dr x dx dip fei, s x ; r, x, <p) x(u; 2 , h, s 2 ; t, x, y?) = 

2tt5(uj 2 - to-i) &T 1 5(A; 2 - fei) 5 S2 _ S1 . (3.25) 

The quantum numbers k > and c<j are continuous, and s is discrete. It takes the values s G Z+A 
or s G 2Z + A, depending on whether the particles are distinguishable or not. This follows again 
from the periodicity condition imposed on the wave function. Taking this into account when 
summing over s, we have the following completeness relation for the eigenstates, 



E 



du k dk x(uj, k, s; r u x 1 , </?i) xO, k, s; r 2 , x 2 , p 2 ) = 

2ti5(t 2 - n) x^5(x 2 - xt) 6\((p 2 - <fi), (3.26) 



where 5\ is again the periodic delta function satisfying the Q3.14| ) or ( p.!5[ ). So, we now have 



a basis of the extended Hilbert space, where the energy and momentum operators are diagonal, 
with eigenvalues 

ft s 

K = ftk, M = ftu, J = fts, S = — . (3.27) 

1 — Mto 

We see that the spectrum of J is the one which is quantized in steps of ft, like the angular 
momentum of the free particles in fll.69p . This is reasonable. On the classical phase space it is 



the charge J which generates a rotation with closed orbits, and a period of 2n. The spectrum of 
S is also discrete, but the eigenvalues additionally depend on the eigenvalues of M. The closer 
the energy gets to the maximal physical energy M max = Mpi/4 = I /AG, the larger the steps are 
between the eigenvalues of S. We shall discuss this spectrum and its physical implications at the 
very end of this section. 

It is now straightforward to impose the mass shell constraint. We can simply repeat all the 
steps from section [I]. The constraint acts on the energy momentum eigenstates as 



£ xd^i k, s; r, x, ip) = (snh 2 (ftk) — A 2 F(ftuj)J xi^i k, s; r, x, (p). (3.28) 
It is again useful to introduce a dispersion relation 

^^^^(^my (3 . 29) 
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It is the quantum version of the classical relation between the momentum K and the energy M, 
which is shown in figure |5|. We conclude that the physical states are those energy momentum 
eigenstates, where the quantum numbers satisfy cu = cu(k). These states are annihilated by the 



constraint fl3.28| ). The most general physical state is a superposition 

V>(t, x, <p) = / k dk i[)(k, s) x(w(k), k , s , T , x , V?), (3-30) 

s J 

where ip(k, s) is again the wave function in momentum space. Let us once again split off the 
dependence on the radial coordinate, writing 

il>(r,x,(p) = Y, I kdke- iu)ik)T e isv *fj{k,s)({k,s]x). (3.31) 

s 

The radial wave function is then given by 



Here we replaced the rescaled hyperbolic functions by the usual ones, to see that only the Planck 
length £ appears as a dimensionful constant. In the limit £ — > 0, we still recover the free particle 
expression ( |1.72| ). 

Position representation 

Before we read off any physical information from the wave functions, we have to transform to 
a representation where the geodesic distance R between the particles is diagonal, and not the 
auxiliary variable X. So, the task is to find the eigenstates of R and then to transform the wave 
functions into this representation. According to ( |2.17| ), we have R 2 = X 2 + (AttGQ) 2 , thus 



R 2 ij(r } x,ip)= (V-(47r£) 2 |U JU) i/>(T,x,<p). (3.33) 
The eigenfunctions of this operator are modified Bessel functions of the second kind, 



'sinh(r/4£) / x 

r](r;x) = - — K jr — — =^ Rrj(r;x) = rr](r;x). (3.34) 



2V2tt 2 £x ire V4tt£ / 
The eigenvalues r are positive, as it should be, and continuous. The eigenstates are orthonormal, 

xdxrj(ri, x) 77(r 2 ; x) = r^ 1 5(r 2 — r{), (3.35) 



and complete 

\drr 1 {r-x 1 )r 1 {r-x 2 )=x-H{x 2 -x 1 ). (3.36) 



The transformation of a wave function from the X-representation to the ^-representation is given 
by 

ip(r,r,(p)— xdxrj(r;x)ip(r,x,(p). (3.37) 
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For simplicity, we use the same symbol for the wave function in both representations. This 
transformation is known as the Kontorovich Lebedev transform $2$j\. The scalar product in the 
^-representation becomes 

(V'llV'a) = / drrdrdv?^i(T,r,</?) V^r,^), (3-38) 

so that the transformed wave function tp(r, r, cp) represents the correct probability amplitude in 
the relative position space of the particles. We can finally write down the most general physical 
state as a wave function in the position representation. All we have to do is to transform the 
radial wave functions (|3.32|) from the X -representation to the ^-representation, thus 



((k,s;r) = J xdx7](r; x) ((k, s; x). (3.39) 

This convolution can be carried out explicitly, and the result can be expressed in terms of a 
hypergeometric function [|25|], 



s T(A)T(A) Jsmh(r/4£) sinh(8vr^) 
L(k,s:r) = -=V"r — -tt ; x 



x sinh A+J? - 1 (4vr£A;) 2 F X (A, A; A + A; - sinh 2 (4vr£fc) ) , (3.40) 



where 

1 / |s| IT \ — 1 / |tS| i7*\ 

A= 2 [ 1+ l-4eu(k) + 4^)> A= 2\ 1+ l-4eu(k)~4^)- ° A1) 

So, we finally have a rather complicated but explicit representation of the physical energy eigen- 
states as wave functions in the relative position space of the particles. The operators for R and (ft, 
representing the distance and the orientation of the particles, are diagonal in this representation. 
It is therefore possible to read off the physical properties of the states directly from the wave 
functions. 



Wave functions 

We shall in the following study the radial wave functions ((k,s;r), representing the physical 
states with eigenvalues 

M = hu(k), K = hk } J — Us, S = ^tt-. (3.42) 

1 - Uu(k) 

We shall compare them to the corresponding radial wave functions C(k, s; r) of the free particle 
system, as given by ( |1.72| ). There we may simply identify the radial coordinate x with r, as in the 
limit G — > the transformation from the X-representation to the ^-representation is trivial. In 
both cases, the physical Hilbert space is spanned by a two parameter family of energy momentum 
eigenstates, and the quantum numbers k and s have the same interpretation. The former is the 
inverse radial wavelength at infinity, representing the momentum of the incoming and outgoing 
particles, and the latter is the inverse angular wavelength, representing the angular momentum. 



48 



Actually, we haven't check this yet for the Kepler system. For the free particles, it is a well 



know feature of the Bessel function ( |1.72| ). For large r it falls off with y/r, oscillating with a 



wavelength of 2ir/k. It is not at all obvious that the function ( P-40[ ) has the same property. But 



we can use the following semi-classical argument. We have seen in ( |2.78| ) that for the classical 
trajectories, the canonically conjugate momentum P of the distance R is equal to ±K for large 
R. Therefore, we expect the eigenfunction of K with eigenvalue Tik to be a superposition of an 
ingoing and an outgoing radial wave, with wavelength 2ix /k for large r. This argument holds for 



both the free particles, where the corresponding classical relation is ( |1.44| ), and for the coupled 
particles. 

An alternative way to see this is as follows. This argument can be made rigorous by taking 
into account the correct operator ordering. But for simplicity let us stick to a semi-classical level. 
Consider the eigenvalue equation for K in the ^-representation. As a phase space function, we 
have 

. Q 2 + S 2 R 2 tnh 2 P + S 2 S 2 - /nA ^ 

Hence, for large R we have P ~ ±K. Since P is canonically conjugate to R, and thus repre- 
sented by the operator —ihd/dr, it follows that an eigenfunction of K with eigenvalue %k has 
a radial wavelength of 2tt /k for large r. At this point, we see what the advantage of the aux- 
iliary variable X was. Of course, we could have used R and P as phase space variables from 
the very beginning. Then, however, we had to solve the eigenvalue equation for (P-43|), with 



the appropriate operator ordering. Obviously, this is much more difficult than solving the eigen- 
value equation for ( |3.21| ), because the hyperbolic functions become a non-local operators in the 
complex r -plane. 

But now, let us come to the actual physical questions. We wanted to find out what the basic 
differences are between the radial wave functions ((k,s; r) for the free particles and those for 
the coupled particles. So far, we know that both represent a scattering state, where the incoming 
and outgoing particles have a radial wavelength of 2n/k, and a fixed angular momentum. In 
other words, the wave functions are almost equal for large r, where large means some orders of 
magnitude above the Planck length i. The interesting question is what happens at small distances. 
To see this, we have plotted some typical wave functions, for different values of k and s. 

In each of the following figures, the value of s is fixed, and that of k varies. The upper part 
of the figure always shows the free particle wave functions, the lower one those for the Kepler 
system, both as a function of the distance r, which is measured in units of the Planck length i. 
To distinguish the wave functions for different momenta k, we use the following rule. The larger 
k is, the smaller are the gaps in the curves. The broken curve with the largest gaps represents the 
smallest momentum, and a solid curve corresponds to the largest momentum in the figure. The 
quantum numbers k are given in units ofl/£, so that k = 1 corresponds to a radial wavelength 
of 2tt£ at infinity, and to a momentum in physical units of 1/G = M P1 . 

The mass parameters are given in units of the Plank mass Mpi. They only have a marginal 
influence on the qualitative behaviour of the wave functions. We therefore keep them fixed in 
the following, and we choose two particles with the same mass m\ = m 2 = 0.02Mpi, thus still 
far below the maximally allowed rest mass. For simplicity, we shall also restrict to states with 
integer s, thus particles with conventional statistics. There are no principally different effects 
occurring for anyons, only the numerical values of some quantities considered below are slightly 
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Figure 6: The wave functions for s = 0. The coupled particles avoid being at the same point 
in space, but otherwise the wave functions are almost equal. 



different. We shall briefly comment on this at the appropriate place. 

As a first example, consider the states with s = in figure In the upper part we see the 
Bessel functions J (kr). They take a maximum at r = 0, and fall off with the square root of 
r for large r, oscillating with a wavelength of 2ix /k. Except for a small region near the origin, 
which is of the order of ten Planck lengths, the wave functions of the Kepler system are almost 
the same. For a vanishing angular momentum, the coupled particles behave almost like the free 
particles. There is not even a phase shift between the ingoing and outgoing particles, so that far 
away from the origin the interaction of the particles cannot even be detected. However, the wave 
function of the coupled particles goes to zero at r = 0, so that the probability to find the two 
particles at the same point in space vanishes. 

Obviously, the coupled particles try to avoid being both in a region of space which is of the 
order of the Planck scale. This becomes even more obvious when we look at the states with 
s = 1 in figure [7| For the free particles, we now have the Bessel functions J\{kr). They are 
zero at r = 0, and the first maximum is approximately at r rs 1/fc. For large r, they still fall off 
with the square root of r, oscillating with a wavelength of 2-k /k. Apart from a different phase 
of the oscillation, this is also the large r behaviour of the wave functions of the Kepler system. 
To see this more clearly, we have plotted the same wave functions for large distances in figure [8|. 
The phase shift indicates that an interaction between the ingoing and outgoing particles has taken 
place. 

The more interesting feature is again the behaviour of the wave function for small r. For the 
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Figure 7: The wave functions for s = 1. The free particles can still get arbitrarily close to 
the origin, whereas the coupled particles alway stay beyond the minimal distance -R m i n - 



free particles, the maximum of the wave function is getting closer to the origin with increasing 
momentum k. For the coupled particles however, the maximum never gets beyond a certain 
minimal distance R m i n - In fact, we already found such a minimal distance at the classical level. 
Remember that the quantum number s specifies the eigenvalue hs of the canonical angular mo- 
mentum J. On the other hand, we have seen that the classical phase space function J is restricted 
by ( p. 12[ ). We used this to show that a Schrodinger quantization, in a representation where X 
and 4> are diagonal, is not possible. 

Turning the argument around, it follows that for a fixed J, there is a minimal X, and thus 
also a minimal R, since R 2 = X 2 + (AnGQ) 2 > X 2 . Consequently, we have the semi-classical 
relation 

2tt 2 £ 

R> R m - m (s) = — \s\- (3-44) 

It follows directly from ( |3.12| ) with J = hs. For a fixed angular momentum, there is thus a 
minimal distance between the particles. It is of the order of the Planck length i if the angular 
momentum is of the order of %, and it increases linearly with the angular momentum. The 
numerical factor in ( p.44[ ) depends on the masses of the particles, defining the number < A < 1 
in ( [2.41D . We see in figure [7| that beyond the point -R m i n the wave function is no longer oscillating, 
but falls off exponentially. In the classically forbidden region, the wave function shows the 
typical tunneling behaviour. 

Another peculiar observation is the following. For large momenta k, the wave function of the 
coupled particles freezes at small distances r. More precisely, for small r and large k the radial 
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Figure 8: For large distances from the origin, the wave functions of the free and the coupled 
system are always the same, except for a phase shift between the ingoing and outgoing wave. 



wave function ((k, s; r) no longer depends on k. The limiting wave function is shown as the 
dotted curve in figure [7]. This is quite strange, because it means that the short distance behaviour 
of the wave function will not change anymore with increasing momentum, whereas for large 
distances we still have an oscillation with a wavelength of 2n/k. One can show that the limiting 



wave function is rj(r; R mm ), as defined in fl3.34h . It is the /^-representation of the eigenstate of 



X with eigenvalue x = -R mm - Hence, for small r and large k we have 

C*(k)(k,s;r,<p)*j»<'x(r;Runn), (3-45) 



with i? m in given by ( p.44[ ). Let us only briefly sketch how this can be proven. If we return 



from the ^-representation to the X-representation, then the radial wave functions are given by 
the Bessel functions (|3.32|). In the limit k — > oo, one can show that these functions provide an 



approximation of the delta function 8(x — R m i n )- Transforming again to the /^-Representation, 
this implies fl3.45| ). But let us not go into any details of the proof here. Numerically it can be 



verified easily that ( |3.45| ) is in fact a good approximation. 

What does it mean physically that the wave function freezes at small r, and that the first max- 
imum never goes beyond the point r = -R mm ? Consider again the scattering process described 
by these wave functions. In case of the free particles, we can make the following statement. 
The larger we make the momentum of the incoming particles, the closer the particles get to each 
other. In other words, to probe small distances we need large momenta. This is a well known 
rule in elementary particle physics or actually quantum physics in general. It is one way to ex- 
press Heisenberg's uncertainty relation. However, this statement is no longer true when gravity is 
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switched on. We can make the momentum as large as we like, we never get beyond the minimal 
distance R mm . Apparently, there is some repulsion that keeps the particles apart. 

To see that this is a real quantum gravity effect, observe that R mm is of the order of the Planck 
length, which involves both G and %. Hence, the effect disappears both if we switch off gravity, as 
then we have the wave functions of the free particles, and it also disappears in the classical limit. 
There is no restriction on the integration constant R for the classical trajectories in section |2|, 
thus the particles can reach any arbitrarily small distance. But what does it mean that the minimal 
distance -R m i n depends on the quantum number s? First of all, one could argue that in order to 
bring the particles closer together we just have to consider states with smaller s. And in fact, we 
have seen in figure that no such minimal distance exists for s = 0, in agreement with ( |3.44[ ). 



But now assume that the particles are two identical fermions with masses m 1 = m 2 = m. 
Then we have s E 1 + 2Z, and consequently the minimal value of \s\ is 1. In this case, we have 

2vr 2 £ 2vr 2 £ 
R > Rmin = —- = — n r , (3.46) 
A cos(47rGm) 

and this is an absolute minimum for the distance, which can never be reached by any wave 
function. This absolute minimum only depends on the mass m of the particles. For m = we 
have -R m in = 2ti 2 £ « 20£. With increasing mass it increases unboundedly. For m — > M P1 /8, 
which is the maximal mass of two identical particles, the cosine goes to zero. For anyons the 
absolute minimum is smaller, because the minimal \s\ lies between zero and one, and for bosons 
this curious effect disappears. 

So, we find the following remarkable physical effect of quantum gravity. If the statistics 
parameter A is different from zero, then the particles are no longer able to approach each other 
closer than some minimal distance, which is of the order of the Planck length. But this is not 
the only feature of the wave functions, from which we can learn something about the short 
scale structure of spacetime. A second effect also results from the relation ( |3.44[ ), and this is 



even independent of the statistics. Consider the wave functions for larger angular momenta, for 
example s = 5 in figure ^. For large r, outside the figure, the free particle wave functions and 
those of the coupled particles are still of the same form, just with a phase shift, as in figure [8|. 

The minimal distance -R m i n is now however five times as large as in figure [7]. Now, suppose 
that we want to localize the particles relative to each other, not at the origin, but at some finite 
distance R = r , and with some angular orientation (p — cp . Clearly, for this purpose we need a 
superposition of states with different quantum numbers k and s. We need large quantum numbers 
s to get a sharp peak in the angular orientation <\>. And we need large quantum numbers k to get a 
sharp peak in the distance R. However, we have just seen that for large angular momenta s, there 
is also a large minimal distance .R mm . And moreover, for large momenta k the wave function 
freezes when the distance is of the order of this minimal distance. 

Roughly speaking, we can say that once we have a sharp angular direction of the particles, 
we do no longer have enough radial states to superpose, in order to get a sharp distance R. Of 
course, this is a somewhat heuristic argument, so let us make it more precise. Consider first the 
following analogy. Suppose that we have a box potential in flat space, which is zero outside 
and positive inside the box. There are then two types of energy eigenstates, those tunneling 
through the box, and those with sufficiently high energy passing the box without tunneling. Here 
we have the situation that actually only the tunneling states exist, because the region r < _R min 
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Figure 9: The wave functions for s = 5. The minimal distance R m i n increases linearly with 
s. As a consequence, it is not possible to localize the particles relative to each other, by 
superposing states with different k and s 



is classically forbidden for all physical states. The physical Hilbert space only consists of the 
tunneling states. 

But this subspace does not provide a complete basis of states in the position space of the 
system. More precisely, it is not possible to superpose the states, so that an arbitrarily sharp 
peak in position space comes out. To make this a little bit more explicit, let us try to define a 
physical state, which describes a localized pair of particles. In case of the free particles, there 
is no problem. In the Schrodinger quantization, we simply have to consider the special wave 
function 

ip(r, <p) = r~ l 5(r - r ) 5\{y - <fo). (3.47) 

It can of course be written a superposition of the energy momentum eigenstates. We can even 
say explicitly what the corresponding wave function ip(k, s) in momentum space is, which has 
to be inserted into Q1.63P - It is 



^(k,s)=e- is ^C(k,s;r ), (3.48) 



because we have the completeness relation ( |1.6ip . The corresponding wave function in the Dirac 



quantization is obtained by inserting the same wave function ip(k, s) in momentum space into the 
formula ( |1.7ip . The resulting wave function tp(r, r, cp) is equal to Q3.47P for r = 0, thus describ- 



ing a localized pair of particle at a given moment of time. As this is not an energy eigenstate, the 
particles will of course be localized only at one moment of time. 
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So far, this is just ordinary quantum mechanics of free particles. But now, let us ask whether 
such a localized state also exists for the coupled particles. Are there physical states describing 
a localized pair of particles, with a wave function of the form ( |3.47| )? We could try to insert the 



same wave function in momentum space ( [3.48D into the general superposition ( B.31D of physical 
states. Doing so, we get 

x/jIt, r, (p) = I k dk e~ iul{k)T ^f-vo) ^ s; rQ ) s; r ). (3.49) 

s J 

Now, what about the functions ((k,s; r) for fixed s? Do they define an orthonormal basis of the 
space of all radial wave function? If this was the case, then we would immediately recover the 
wave function Q3.47D for r = 0, thus we have a localized pair of particles at that moment of time. 



Unfortunately, however, the functions ((k, s; r) do not provide an orthonormal basis. This is not 
obvious from the representation fl3.40j ), but it can be seen by going back to the X -representation 
(J3.32p. The Bessel functions provide an orthonormal basis if and only if the index is fixed. For 



the free particles this is the case, since in ( |1.72[ ) the index is \s\. 



In 03.4OD , however, the index is |s|/(l — A£u(k)), and thus not only the argument but also 
the index of the Bessel function depends on k. The functions ((k; s, r) for fixed s in ( |3.49j ) are 



not orthonormal, and therefore the resulting wave function for r = is not of the form ( |3.47| ), 
describing a localized state. Now, one could argue that there might be another superposition of 
physical states describing a localized pair of particles. However, if there were such localized 
states, then we could use them as a basis of the physical Hilbert space, and this would be exactly 
the basis that we defined in the beginning of this section, when we tried to set up a Schrodinger 
quantization. We saw, however, that due to the global structure of the classical phase space such 
a basis does not exist. 

So, what is the conclusion? Obviously, the spacetime in which the quantized particles are 
moving has some kind of foamy structure. Maybe it is possible to analyze this structure more 
explicitly, by having a closer look at the wave functions ((k, s; r). For a single particle system, 
such an analysis has been carried out in some detail in [|T3"|], finding a semi-discrete structure of 
spacetime and an explicit uncertainty relation, which forbids the particle to be localized. Here 
things are more involved and therefore it is more difficult to make the idea of a foamy spacetime 
more precise. But at least we can see that it is impossible to localize the particles at a point in 
space, and this is already a nice result. It gives us a hint to what kind of influence a real theory 
of quantum gravity might have on the local structure of spacetime in general. 



The quantized cone 

Finally, let us consider yet another strange result, which has actually not very much to do with 
the point particles. There is quantized conical geometry of the spacetime at spatial infinity. Let 
us once again look at the spinning cone (|2.7D. It defines the asymptotic structure of the spacetime 
at infinity. We saw that it has the following global geometric features. There is a deficit angle of 
8nGM, and a time offset of 8ttGS. Now, consider the energy momentum eigenstates spanning 
the physical Hilbert space. They are the eigenstates of M and S, with eigenvalues 

Tis 1 
M = hu } S = - — , 0<lo<— , s£Z. (3.50) 
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For simplicity, we assumed that the particles are massless, distinguishable, and that the statistics 
parameter is A = 0. There is nothing essentially different in a more general case, except that the 
formulas below are slightly different. In the context of quantized general relativity, the energy 
momentum eigenstates can also be regarded as states with a sharp geometry of the spacetime at 
infinity, but with the locations of the particles in space smeared out. However, the geometry of 
the spinning cone is thereby not arbitrary. There is an obvious relation between the eigenvalues 
of M and S, namely 

(1 -AGM)S e Zh. (3.51) 

As this results from the quantization of the canonical angular momentum J, we expect that this 
is a general result, which does not depend, for example, on the number of particles present. In 
fact, one can show that for a general multi particle model there is the same relation between the 
canonical angular momentum J and the parameters of the spinning cone M and S, and the same 
is expected for any kind of matter included [ ]T0| ] . 

Thus, we conclude that the geometry of the conical infinity itself is quantized. Just like in 
a quantized atom, where not every classically allowed combination of energy and angular mo- 
mentum can be realized, the conical geometry of the spacetime at infinity is restricted by this 
quantum condition. It becomes a very simple relation when it is written in geometric units, that 
is lengths and angles. We define the total angle of the spinning cone to be a = 2ir — 8ttGM. It 
is the circumference of a unit circle around the central axis. And we define the time offset to be 
r = 8nGS. Then the quantization condition becomes 

arG(47r) 2 Zl (3.52) 

Only those conical spacetimes can be realized where the product of the total angle and the time 
offset is an integer multiple of (47r) 2 £. Note that a is dimensionless and r is a length, or time, 
so that necessarily the Planck length shows up. If a is close to 2%, then r is quantized is units 
of 8n£, which is a kind of time quantization. If we assume that this result is independent of the 
kind of matter being present inside the universe, then what we have found is a real quantization 
of geometry. 
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Appendix 

Here we summarize some facts about the spinor representation of the three dimensional Lorentz 
algebra sl(2) of traceless 2x2 matrices, and the associated Lie group SL(2). A more compre- 
hensive collection of formulas, using the same notation, can be found in []I3[], As a vector space, 
fil(2) is isometric to three dimensional Minkowski space. An orthonormal basis is given by the 
gamma matrices 

7o=(_° 1 J), 7i= (J I), 12= (J \). (A.1) 
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They satisfy the algebra 

lalb = Vab 1 - ^abc Y (A.2) 

where a, b = 0, 1, 2, the metric rj a b has signature (—,+,+), and for the Levi Civita symbol e abc 
we have e 012 = 1. Expanding a generic matrix in terms of these gamma matrices, we obtain an 
isomorphism of sl(2) and Minkowski space, 

v = v a la v a = \Tv{v 1 a ). (A.3) 

Some useful relations are that the scalar product of two vectors is equal to the trace norm of the 
corresponding matrices, and the vector product is essentially given by the matrix commutator, 

\Ti{vw) = v a w a , \[v, w] = -e abc v a w hlc . (A.4) 

Sometimes it is useful to introduce cylindrical coordinates in Minkowski space, writing 

v = rjo + p7(v?), (A.5) 

where r and p > are real, and (p is an angular direction. The vector 7(93) and its deriva- 
tive j'(p) form a rotated set of spacelike unit vectors, pointing into the direction of p and the 
orthogonal direction, 

7(<^) = cosy? 71 + sin (/? 72, 7' ((p) — cos ^72 — sin ip 71. (A. 6) 

Useful relations are 

7o7(y?) = TV), 7o7 / (y) = -7(v), ( A - 7 ) 

and 

7(y?i) 7(^2) = cos(v?i - ip 2 ) 1 + sin(v?i - ip 2 ) 7o- (A.8) 

The Lie group SL(2) consists of matrices u with unit determinant. The group acts on the algebra 
in the adjoint representation, so that 

v 1— > u l vu. (A.9) 

This provides a proper Lorentz rotation of the vector v. The conjugacy classes of the algebra 
are characterized by the invariant length |Tr(t> 2 ). For timelike vectors v with ^Ti(v 2 ) < 0, we 
distinguish between positive timelike vectors with v° = ^vj° > 0, and negative timelike vectors 
with v = |«7° < 0. And similar for lightlike vectors. 

A special group element, which represents a clockwise rotations about the 7 -axis by an angle 
a, is 

u = e a ^ ,/2 = cos(a/2) 1 + sin(a/2) 7o . (A. 10) 

This implies 

70 ^ -1/^70 u = 70, j(p) ^ u~ 1 j(p)u = j(p - a). (A.ll) 
A boost is specified by a rapidity x > and a direction (3, 

u = e x ^ = cosh X 1 + sinh x 7 (/?) . (A. 12) 
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A generic element u 6 SL(2) can be expanded in terms of the unit and the gamma matrices, 
defining a scalar u and a vector p 6 sl(2), 

u = Ml+/7 a p = p a la . (A.13) 

The scalar w is basically the trace of u, and the vector p is called the projection of u. The 
determinant condition implies that 

u 2 =p a p a + l. (A. 14) 

According to the property of the vector p, we distinguish between timelike, lightlike and space- 
like group elements u. For a timelike element we have — 1 < u < 1, and it represents a rotation 
about some timelike axis, which is then specified by the vector p. The angle of rotation a is 
given by u = cos(a/2), and the direction is given by the sign of p°. 
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